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(HHC) flight test program showed significant reductions in main rotor shaft torque 
and engine power in hover and forward flight [1]. The unsteady aerodynamics with 
the higher harmonic control svstem application, including wake etfects. were con- 
sidered to studv whether such power reductions are feasible. 

An airtoil oscillating in pure plunge can achieve propulsive torce (" Katzmavr 
effect"), as in the case of birds’ wings flapping in flight. Here it will be shown 
that this effect can be enhanced in the presence of layers of shed vorticity with 
the proper phasing. In addition, while it is known that an airfoil oscillating in 
pitch, will typically produce drag at most values of reduced frequency, it is found 
that the presence of another layer of shed vorticity of the proper phase, can 
reduce the drag on the pitching airfoil depending upon wake spacing, reduced fre- 
quency, and phase. Under some conditions the added layer or layers of shed vor- 
ticity will even result in propulsive force acting on the pitching airfoil similar to 


the "Katzmayr" effect for the plunging case. 

It was found, for the OH-6A helicopter, that the measured reductions in 
main rotor shaft torque and engine power are feasible when evaluated with respect 
to the "Katzmayr effect" and the additional drag reduction or propulsive force 
obtained due to pitch and plunge oscillations with the effect of adjacent wake 


layers of shed vorticity. 


sea 


TABLE OF CONTENTS 


1. INTRODUCTION 


A. GENERAL هس‎ 200 | 
B. ROTARY WING UNSTEADY AERODYNAMICS .........›.......›››› D 
C. ADDITION OF WAKE EFFECT COTES 5 
D. HIGHER HARMONIC CONTROL ( HHC ) .......... cernentes 6 
E. OH-6A TEST PROGRAM DATA مس ٗػٗصك>ٗپ”صأأى٭پاػػممممسب‎ 7 
F.SCOPE ............. ne s 12 
II. BACKGROUND RL nun ee TEE AS 15 
III. DRAG REDUCTION DUE TO WING FLAPPING 
(KATZMAYR EFFECT ) 
A. INTRODUCTION رر سے سے سس سس سو‎ 31 


IV. 


. IHE. PHYSICS OF FLAPPING AIRFOILS IN 


PON ٢٢۱٢ ٧٧٢ AAA A E Ena ce nhi dee 32 
C. PANEL CODE FOR STEADY INCOMPRESSIBLE 

IS DON een ne ae 34 
DANA AO — . aoreegeevesgesceteCeascteetuee use ee ete oreet 42 
الا‎ ODOIRSENSSENNONOY SIS iie ee 1 SSS 46 
F. GARRICK'S PROPULSIVE FORCE ANALYSIS .................... 47 
ሸባብ እ11.17)ካ (71. [፡1..31 11.3 5e aiieces te القن ما جم م لامع لع عام ولمع ما الال‎ 48 
ll. ESTIMATED POWER REDUCTION DUE TO IHC 


O SADA AD A E N N tes access SD 22 


ANALYSIS OF WAKE INTERFERENCE EFFECTS 


X. 


B. 


C. 


D. 


60 یی 0000 10 000080800010010 0 0 E A‏ 
61 الما لل لماه ...2 UNS DENDIZIIN O-KOTTZBPANEL TEODE‏ 
63 مہو مھا ہم مھ سای سس دو د وو واو ۳۳۱۱۲۷۰۰۰۰٠٦٠۹٦۰‏ 
EVALUATION OF TWO-FOIL PANEL CODE ........................ 68‏ 


E. ANALYSIS OF OSCILLATORY BLADE 


INTERFERENCE EFEFEGCTS. 5.5. ee AAEE 74 


V. PROPULSIVE EFFICIENCY AND FLUTTER ANALYSIS 


106 د مر د سای .... A. INTRODUCTION‏ 
107 ادن ند د مس ا ع —.............. B. PROPULSIVE EEEICIENCY‏ 
C. SINGLE-DEGREE-OF-FREEDOM AIRFOIL FLUTTER .......... 115‏ 


VI. RESULTS AND CONCLUSIONS 


A. DISCUSSION OF RESULTS ............... 7ء مم‎ 120 
B. CONCLUSIONS ...uessesnteic venum UU LT ک‎ 132 
C. RECOMMENDATIONS .......... ل ا‎ ٦ 134 
List of References 00101011 eee 135 
Initial. Distribution9list ...... سس سس‎ mE 142 


vi 


Table (5.1) 


Table (3.2) 


Table (3.5) 


Table (4.1) 


Table (6.1) 


LIST OF TABLES 


Page 
The OH-6A helicopter rotor parameters 53 
Cd ys እቢ tor NACA 0015 plingine with hs — |" we... 35 
Cdavepcersussneduccdeimequencv o Kp al NOG 1 os 59 
Conectonerac ROS OPEM TDI DG EE ce occ መ ው 72 


Aerodynamic coefficients vs non-dimensional time (from 


time step 113 to 159) for the case of wake phasing , m =]75 


degrees. wake spacing h=2, reduced frequencv k =0.0617 


29[ مات سس ل2 11603 01 and DIMERS‏ 


Vil 


LIST OF ILLUSTRATIONS 


Fina. 1.1) L/Lo phase diagram ( Imaginary vs Real ) 

Fiuk Power reduction results from higher harmonic control 
open loop test program (NASA/Army/Hughes - Sept. 1982) 

Figure ( 1.3 ) Comparison of 4/rev vertical vibration at the pilot seat versus 
airspeed with HHC "on" and "oft". 

Figure ( 1.4) Comparison of 4/rev. lateral vibration at the pilot seat versus 
airspeed with HHC "on" and "off". 

Figure (2.1) Loewy’s schematic elements of unsteady rotor flow field 

Figure (2.2) Schematic representation of unsteady rotor flow field. 

Figure (2.3) Two-dimensional models of unsteady rotor flow. 

Figure (3.1) Propulsive force on plunging airfoil 

Figure (5.2) Panel Method Geometry 

Figure (3.3) Relationship between geometrical quantities 

Figure (3.4) Panel method representation for unsteady flow 

Figure (3.5) Comparison of panel code results with Garrick ( k=0.1 ) 

Figure (3.6) Comparison of panel code results with Garrick ( K 71 ) 

Figure (3.7) Comparison of panel code results with Garnck ( k=2 ) 


Figure (3.8) Blade deflection vs radial station at 2P 


viii 


Breure (53.10) 


ute ( 5.11) 


elie ( 3.12) 


Figure ( 4.1 ) 


Figure (4.2 a) 


Figure ( 4 2 b) 


Figure (4.3 a) 


Figure (4.3 b) 


Figure (4.4 ) 


Figure (4.5 a) 


Figure (4.5 b) 


Figure (4.6 ) 


Figure (4.7) 


Blade deflection vs radial station at 4P 

Blade deflection vs radial station at 5P. 

Cdav vs kp ( plunging amplitudel" ). 

Aerodynamic unsteadv model for a multi-bladed rotor. 
CLES 16 ne Codes) 

CI/Cl, vs h/c ( From Ref. [44] ). 

Local velocity distribution; V,/V vs X/C ( Panel Code ). 
Local velocity distribution, V/V vs X/C (after Wegley [46]). 
Correction factor, B vs the chord-gap ratio h/c. 

Pressure distibution vs X/C (Panel Code ). 

Pressure distnbution vs X/C (Panel Code ) 

Wake pattern (airfoil #2 plunging at high reduced frequency 


and close to airfoil #1). 


Wake interference between two plunging airfoils in closer 


proximity. 


Figure (4.8) 
Figure (4.9) 
Figure (4.10) 
Figure (4.11) 
Figure (4.12) 
Figure (4.13) 
Figure (4.14) 
TELS) 
Figure (4.16) 
Figure (4.17) 
Figure (4.18) 
Figure (4.19) 
Figure (4.20) 
Figure (4.21) 
Figure (4.22) 
Figure (4233) 
Figure (4.24) 
Figure (4.25) 
Figure (4.26) 
Figure (4.27) 
Figure (4.28) 
Figure (4.29) 
Figure (4.30) 
Figure (4.31) 


CI vs T* ( m= 0, h= 200). 
Cdvs T* ( m=0, h=200). 
Wake pattern ( m=0,h= 2) 

Cl vse 00 

Cd vs T* ( m=0, h=2 ). 

Wake pattern (m=.20833, h=2). 
Cl vs T* ( m= 20833, h=2). 
Cd vs T Em ۰٦ 
Wake pattern im=025.h=2) 

Glevs: T* sem ٢٢ በከ ን.) 

Cd vs T* (m2025 , h22). 

Wake pattern (m=0.5 , h=2). 

Cl vs T* (m=0.5 , h=2). 

Cd vs T* (m=0.5, h=2). 

Average propulsive force coefficient versus wake phasing. 
Wake pattern ( pitch, m=0, h=2 ). 

Cl vs T* ( pitch, m=0, h=2). 
Cd vs T* ( pitch, m=0, h=2). 
Cl vs T* ( pitch, m=0.20833, h=2). 
Cd vs T* ( pitch, m=0.20833, h=2). 
Cl vs T* ( pitch, m=0.25, h =2). 
Cd vs T* (pitch,m=0.25 , h=2). 
Wake pattern ( pitch, m=0.5, h=2 ). 


Cl vs T* ( pitch, m=0.5,h=2). 


Figure (4.32) 
Figure (4.33) 
Figure (4.34) 
Figure (5.1) 
Figure (5.2) 
Figure (5.3) 


Figure (5.4) 


Figure (5.5) 


Figure (5.6) 


Figure (5.7) 


Figure (5.8) 


Figure (6.1) 


Cd vs TweDrich, m-05, hz2) 

Average drag vs wake phasing in case of pitch. 

Average drag vs reduced frequencv k. 

Propulsive efficienev 5 versus l/k for a plunging airfoil. 
Propulsive efficiency 9 versus l/K for an airfoil pitching 
about the leading edge. 

Propulsive efficiency 1 versus I/k for an airtoil pitching 
about the quarter-chord point andmid-chord point. 

Propulsive efficiency nn versus l/k tor an airfoil 

pitching. about the three - quarter- chord point. 

An airfoil restrained to rotate about its leading edgein two- 


dimensional flow. 


Variation with reduced frequency k of the real and imaginary parts 


of the non-dimensional aerodynamic moment my due to pitching 


of an airfoil about its leading edge in an incompressible flow. 


Comparison of Cm (Imag.) vs. k, ( NACAO0007 oscillating in 


pitch about the leading edge for a, =1° and N=200 panels ). 


Imag. (mv) vs k (for an airfoil pitching at ea=-0.5, 0.0.0.5, and 1.0) 


Drag coefficient vs angle of attack. 


xi 


NOMENCLATURE 


angle of attack 


pitching amplitude 


alr densitv 


rotor angular speed 


frequency of oscillation 


natural frequency of oscillation in plunge 


natural frequency of oscillation in pitch 


coordinate of axis of rotation 


airfoil half-chord 


airfoil chord 


airfoil lift coefficient 


reference airfoil lift coefficient 


second airfoil lift coefficient 


airfoil drag coefficient 


xii 


Lo, Ln 


Mp, M 


A 


average drag coefficient 


reterence airfoil drag coetficient 


second airfoil drag coefficient 


moment coetficient 


propulsive torce coetficient due to plunge 


propulsive force coetficient due to pitch 


average drag force 


non-dimensional wake spacing 


non-dimensional plunging amplitude 


reduced frequency (based on airfoil half-chord ) 


reduced frequency ( used bv panel code and based on airfoil chord ) 


lift force 


aerodynamic coefficients used for Theodorsen analvsis 


moment 


aerodynamic coefficients used for Theodorsen analysis 


xiii 


m 


Pxh 


ps 


UN 


x 0 


ratio of frequencv of oscillation to rotor speed 


number ot blades 


propulsive force due to plunge 


power 


rotor radius 


local radius 


blade area 


non-dimensional time 


free stream velocity 


coordinate of axis of rotation 


xiv 


Acknowledgments 


This research was conducted at the Department of Aeronautics and Astronautics. 
Naval Postgraduate School, Monterey, California, under supervision of Prof. E. 
R. Wood. Chairman of my doctoral committee, Prof. M. F. Platzer, Prof. G. V. 
Hobson, Prot. D. Danielson (Department of Mathematics), and Prot. J. V. Sanders 
(Department of Physics). 

| would like to express my sincere appreciation. to. Prof. E. Roberts Wood 
who provided me with this great opportunity to work in this area and for his 
continuous support. ideas, guidance, knowledge and encouragement. Special thanks 
and deep gratitude are also due to Prot. Max F. Platzer for his continuous support. 
guidance, knowledge and valuable hours of council. Also I would like to express 
my gratitude and appreciation to Prof. Garth Hobson for his help and support. 
guidance, encouragement, and improving my knowledge in classes and in conducting 
this research. 

Deep appreciation goes to Prof. Donald Danielson for his sincere guidance 
and the great pleasure to discuss and understand the physics behind the mathematical 
problems during the hours of his special and useful Applied math. class. Special 
thanks to Prof. James Sanders and Prof. Allan Coppens for their great help, guidance 
and valuable knowledge they provided me in the Department of Physics. 

Special appreciation is due my mother, my brother and my sister. Without moral 
support of my wife, Neven, my lovely daughter, Yasmin, and my nice son Aboda, 


this work could not have been carried out. To them I dedicate this thesis. 


Summary 


Open loop performance data trom the NASA-Army OH-6A higher harmonic 
control (HHC) flight test program. Wood et al. [1], show significant reductions 
in main rotor shaft torque and engine power in the airspeed regime trom hover 
to 100 knots. Depending upon HHC controller phase and helicopter airspeed, 
reductions in power were recorded as large as 2096, with reductions of the order 


of 10% being typical. 


Higher harmonic control 1s an active helicopter vibration control concept 
which alters the aerodynamic loads on the rotor blades such that the blade re- 
sponse is reduced. This in turn reduces the vibratory forces and moments acting 
at the hub, which cause airframe vibration. Basically, HHC is an electronic, 
computer-controlled active vibration suppression system which senses and cancels 
vibrations in a helicopter by N/revolution feathering or pitch motion of the rotor 
blades. N being the number of blades'. Very substantial vibration reduction was 
achieved under a NASA-Army sponsored program using a modified OH-6A he- 
licopter. Higher harmonic blade pitch control was implemented by superimposing 


4/rev. (32 Hz) swashplate motion upon basic cyclic and collective control inputs. 


l E. Robens Wood, R. W. Powers. J. H. Cline. and C. E. Hammond. " On Developing. and 


Flight Testing a Higher Harmonic Control Svstem." Journal. of A.H.S Januarv. 1985. 
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During the program. in addition to reducing vibration levels, the HHC system 
showed the potential for decreased helicopter power requirements. Recent studv 
of the OH-6A flight test data indicates that the mechanism bv which the power 
reduction was achieved is related to the unsteady aerodynamics associated with 
HHC. This is supported by the fact that higher harmonic control by the very 
nature of the method by which it is implemented through oscillating an airfoil 
in pitch with resulting plunge motion, requires definition of unsteady aerodynam- 


Ics to adequately model its physics! 


The fundamental closed form solutions of Theodorsen [3], and Loewy [4] 
provide the basis for theoretical work in this area. The closed form theory 
shows rapid changes in the lift deficiency function with changes in reduced fre- 
quency, wake spacing and frequency ratio. In the past. emphasis in the study 
of unsteady aerodynamics has tended to focus on flutter instability and the effect 


of unsteady aerodvnamics on generating lift and torsional loads. 


In this research, we are interested in performance. and emphasis will be on 
the effect of unsteady aerodynamics on the drag of the airfoil. When the drag 
force is reversed, acting to propel the airfoil forward. it is generally referred to 
as either negative drag or propulsive force. The classic reference on this subject 


is that by Garrick [5]. 


! E. R. Wood. Max F. Platzer. Ahmed Abourahma. Mark A. Couch " On the Unsteadv 
Aerodynamics of Higher Harmonıc Control." Paper No. C 17. Nineteenth European 


Rotorcraft Forum. Cernobbio (Como), Italv, September 1993. 


While Garrick’s [5] work shows that an airfoil oscillating in pitch will 
typically produce drag in the lower reduced frequency range, k. it 1s found that 
the presence of another laver of shed vorticity of the proper phase can reduce 
the drag on pitching airfoils and depending upon wake spacing, reduced frequen- 
cy, and phase. may even enhance the propulsive torce acting on the airfoil similar 


to the "Katzmayr" effect for plunging airfoils. 


The specific mechanisms, that have been used in this research to discuss 
the phenomenon of OH-6A helicopter rotor power reduction. are: (1) the unsteadv 
..odynamics associated with a plunging airfoil and how that produces a propul- 
sive force or decrease in rotor power; (2) the influence of a layer of wake 
vorticity on this propulsive force, which represents the influence of the unsteady 
aerodynamics of one rotor blade on adjacent rotor blades, (3) the unsteady aero- 
dynamics associated with an airfoil oscillating in pure pitch and how this pro- 
duces in general a drag force on the airfoil unless accompanied by layers of shed 
wake vorticity ; and (4) the influence of the phasing of a layer of shed vorticity 
on this drag force, and how. depending upon specific phasing, the effect on the 


airfoil will be either drag or propulsive force. 


The analysis of the above effects was based upon: (1) the classical flat 
plate unsteady aerodynamic theory for an incompressible fluid of Theodorsen [3], 
its extension to include chordwise forces (drag) by Garrick [5], the rotor wake 
analysis of Loewy [4], (2) the panel method developed by Platzer et al. [6], 
that can be applied for unsteady incompressible flow past airfoils or airfoil com- 


binations of arbitrary geometry and amplitudes of different kinds of motions. 


xviii 


in the first approach. Theodorsen’s well known lift deticiency tunction C(k) 
has been modified by Loewy for the case of a hovering rotor. Airfoil drag or 
propulsive torce was determined by the methods of Garrick [2]. In the second 
approach. the panel code allowed a systematic study of the lift and drag produced 
by oscillating airfoils and airfoil combinations in an incompressible inviscid 


flow. 


The panel code results were checked against different references. The results 
showed good agreement with these references. For the parameters ot the OH- 
6A rotor in hover, the main conclusion was that the reductions in drag and 
propulsive forces generated on an airfoil, oscillating in pitch or plunge, due to 
wake phasing, are sufficient to explain the reductions in power recorded on the 


OH-6A during the NASA-Army flight test program [1] 
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I. INTRODUCTION 
A. GENERAL 


The determination of the unsteady aerodynamics associated with a helicopter rotor has 
presented a challenge to the helicopter analyst since the conception of the helicopter. The 
etfect of blade oscillations on resulting airloads has prompted new interest in oscillatory 
blade aerodynamics in connection with possible performance gains. Rotary wing unsteady 
aerodynamics is considerably more complex than fixed wing aerodynamics. 

The simplest case is that of hover, but even here the problem is complicated by the 
influence of the shed helicopter wake from the rotor. In forward flight, the problem is 
further complicated by the fact that the velocity of a local blade element ranges at the 
blade tip from transonic on the advancing side to low subsonic on the retreating side, and 
at the blade root from low subsonic on the advancing side to reversed flow on the retreat- 
ing side. 

For this reason there is a great need for advances in the present theorv :n order to 
explain modern helicopter problems. yet, the number of advances in rotorcraft unsteady 
aerodynamic theory are very limited, so we need to review current unsteady aerodynamic 
theory for rotorcraft and to advance the theory in one of the areas where improvement is 
badly needed. 

The goal of this thesis is to explain the measured power reductions in the OH-6A 
higher harmonic control (HHC) flight test program conducted by E.R. Wood et al [1]. 
This research will consider the unsteady aerodynamics with HHC system application, 


including wake effects. in order to determine whether such power reductions are feasible. 


1 


B. ROTARY WING UNSTEADY AERODYNAMICS 
In this section. some elements of unsteady airfoil aerodynamics will be discussed. 
Consider an airfoii undergoing sinusoidal pitch and plunge motion. As the airfoil performs 
an oscillation. vortices are shed into the medium with a circulation strength equal in mag- 
nitude to the increase in circulation about the airfoil, but opposite in direction. These dis- 
turbances are stored in the tluid because the shed vorticity convects downstream at the 
local flow field velocity. 
The counter-rotating vortices induce a sinusoidal flow field which further changes the 
. The result of this complex flow field is a time difference or delay between the 
airfoils motion and the induced aerodynamic forces. This delay is known as the phase 
lag, w [34]. To simplifv calculations for this type of motion, it is common to describe 
the airfoil position and the associated aerodynamic forces by complex variables. For pure 
plunge oscillations the vertical motion of the airfoil is described by the real part of the 
following equation: 
ከ0) = ኪ, ll E 
where h, is a complex number. and © is the frequency of oscillation. Similarly, the lift 1s 


a oD yee 


where L, is the quasi-steady lift given by the expression : 


L5 pU^S CAVA. A دل‎ (1.3) 


This is termed quasi-steady because the angle of attack is represented by h/U. 
The values of r, and w represent the magnitude and the phase respectively, of the 
true instantaneous lift relative to the quasi-steady lift. The variables r, and y in 
general depend on the reduced frequency k, the Mach number M. and the Reynolds 
number. For inviscid. incompressible fluid the values of r, and w will be only a 
function of k [54]. A complete solution for the oscillating tlat-plate airfoil in an 
incompressible inviscid flow has been obtained by  Kussner and Theodorsen and 


reproduced from Reference [35] in Figure (1.1 ). 





Figure (1.1) L/L, Phase Diagram ( Imaginary vs Real). 


C. ADDITION OF WAKE EFFECT 

For the helicopter in hover. the rotation of the wing introduces a number of features 
that require special attention. notably the returning vortex wake. the time varying free 
stream. and radial flow. and a fundamentally transcendental geometry that requires either 
approximate or numerical solutions. The lift on a wing is due to its bound circulation. 
Conservation of vorticity in the flow requires that there be trailed and shed vorticity in 
the wake behind the wing. 

The spanwise variation ot bound circulation results in trailed vorticity parallel to the 
free stream direction. Time variation of the bound circulation leads to shed vorticity par- 
allel to the wing span. The wake is composed of sheets of vorticity convected downstream 
from the trailing edge by the free stream velocity. The wake of a rotor in hover or in ver- 
tical flight consists of helical vortex sheets below the disk, one from each blade. Unsteady 
motion of the rotor blade will produce a shed vorticity in the wake spirals. With low disk 
loading (at low collective pitch settings ), the wake remains near the rotor disk and there- 
fore passes close to the following blades. Thus the wake vorticity is not convected down- 
stream of the airfoil as with a fixed wing, and the shed vorticity sheets below the rotor 
disk must be accounted for to correctly estimate the unsteady loads. 

For high inflow (high collective pitch settings ) or in high-speed forward flight, the 
wake will be convected away from the blades. The returning shed wake influence is pri- 
marily a concem of hover and vertical flight, and forward flight at airspeeds less than 
100 knots. Assuming high aspect ratio of the blade. lifting line theory requires a knowl- 
edge of the loads on the blade section, and the returning shed wake of the rotor must be 


incorporated into the two-dimensional unsteadv airfoil theory. 


The wake far trom the blade section will have little influence. so the emphasis is on 
modeling the wake near the blade. which for low inflow consists of vortex sheets that are 
nearlv planar surfaces parallel to the disk plane. Based on these considerations, a two- 
dimensional model for the unsteady aerodynamics of the rotor can be constructed. Loewy 
in 1957, developed a two-dimensional model for the unsteady aerodynamics of the blade 


in a hovering rotor, including the effect of the helical shed wake. 


D. HIGHER HARMONIC CONTROL (HHC) 

The primary source ot helicopter vibration is the higher harmonic blade loads gener- 
ated by the aerodynamic environment at the rotor disk. For an N-bladed rotor, the oscil- 
latory (N-1)P, NP, and (N+1)P blade loads are transmitted as exciting hub forces and 
moments to the airframe at a dominant frequency of N per rev (NP). A significant reduc- 
tion in vibration level has been achieved by industry in the last twenty five years. The 
desired goal was to achieve 0.02 g and this goal cannot be met without a quantum 
advance in vibration control technology. The vibration levels first specified by the U.S 
Army for the AAH/UTTAS procurement in the early 1970s were later revised upward in 
the mid 1970s to better reflect realistic design goals consistent with the state-of-the-art in 
helicopter vibration control [1]. 

A higher harmonic control system drives the blades in pitch at the (N-1)P, NP, and 
(N+1)P harmonics of rotor rotational speed, generating new unsteady airloads, which in 
combination with the new oscillatory inertial loads, cancel the harmonics of blade loads 
that cause airframe vibration. Therefore the vibrations are suppressed at the source. 
Although there are several ways to implement HHC on a rotor system, the approach 


generally followed to date has been by blade root feathering using swashplate oscillations . 


By means of electro-hvdraulic servo-actuators. the swashplate is excited in the 
collective, longitudinal cyclic. and lateral cyclic modes at NP resulting in blade pitch 
oscillations at three distinct frequencies of (N-1)P, NP. and (N+1)P in the rotating frame. 
In the system considered. developed for the four-bladed OH-6A helicopter. higher har- 
monic blade feathering for vibration reduction is achieved by superimposing 4/rev swash- 
plate motion upon basic collective and cyclic flight control inputs. 

Perturbing the stationary swashplate at 4/rev both vertically and in pitch and roll 
results in third. fourth. and fifth harmonic blade feathering in the rotating svstem. Fourth 
harmonic blade feathering is achieved by oscillating the swashplate vertically about its 
collective position, while third and fifth harmonic blade feathering in the rotating system 
results from 4/rev tilting of the stationary swashplate in pitch and roll about its cyclic tilt 
position [2]. Results of experimental efforts by Hammond [52] showed that successful 
suppression of vibration can be achieved by oscillating the blades at relatively small 


angles. 


E. OH-6A (NASA/ARMY/HUGHES) TEST PROGRAM DATA 

The goal of this research has been a systematic study of the OH-6A helicopter 
HHC flight test program. in order to understand and explain the significant reduction in 
main rotor shaft torque and engine power that was measured in flight. Shown in Figure 
(1.2) are the power reduction results from the HHC Open Loop Flight Test Program 
(NASA/ Army/ Hughes - Sept. 1982). 
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l'igure ( 1.2 ): Power Reduction. Results from Higher Harmonic Control (LEC) 


Open Loop Test Program (NASA/Army/Highes - Sept. 1982) 


Figures (1.3) and (1.4) (From Ref. [1]) show the very substantial vibration. reduction 
achieved under a NASA-Armv sponsored program using a modified OH-6A helicopter. Higher 
harmonic blade pitch control was achieved by superimposing 4/rev. (32 Hz) swashplate 
motion upon basic cyclic and collective control inputs. The aircraft was flown from zero 


airspeed to 100 knots with the LHC system operated both open loop (manually) and closed 


loop (computer controlled). 
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Vigure ( 1.3 ): Comparison of 4/rev Vertical Vibration at the Pilot 


Seat Versus Airspeed with 11116 "on" and "off". 
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Figure (1.4): Comparison of 4/rev. Lateral Vibration at the Pilot Seat 


Versus Airspeed with HHC "on" and "off". 
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The signiticance of power reduction means that: (1) the helicopter pavload can be 
increased. (2) the maneuverability can be markedly increased. (3) marginal operating condi- 
tions, such as at sea in high sea states can be better met. and (4) the helicopter operating 
costs can be substantially reduced. Increased interest in HHC and its potential benefits for 
the helicopter now require that we try touse the proper tool to understand this phenome- 
non. It would appear that we presently have the computational tools, namely the panel 
method for unsteady inviscid, incompressible flow past airfoils or airfoil combinations of 
arbitrary geometry, developed by Platzer et al [6], and the single wake ciosed form anal- 


ysis developed from Loewy’s theory by Wood and Couch [40]. 
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I. SCOPE 


Chapter II contains a quick review of the most important research starting with the 
work of Knoller[10] in 1909 and Betz [11] in 1912. Thev explained the bird's ability to 
generate forward thrust by means of wing flapping. In 1922. R. Katzmayvr, experimentally 
verified the work of Knoller and Betz. Birnbaum, von Karman. Theodorsen. and Garrick in 
the years 1924 through 1936, were the pioneers of flapping wing aerodvnamics. 

The fundamental closed form solutions by Theodorsen [3], Loewv [4], Garrick [5], pro- 
vide our basis for theoretical work in this field. While Loewy's work on wake-induced 
flutter helps to explain the phenomenon. it also points up difficulties to be overcome. That 
is the closed form theory shows rapid changes in the lift deficiency function with changes 
in reduced frequency, k; wake spacing, h; and frequency ratio, m. 

Garrick further developed the work of Theodorsen by deriving the horizontal force for- 
mulas for airfoils or airfoil-aileron combinations oscillating in any of the three degrees of 
freedom: vertical flapping, torsional oscillation about a fixed axis, and angular oscillation 
of the aileron about a hinge. Also presented in Chapter II. Loewy's two-dimensional model 
for the unsteady aerodynamics of the blade of a hovering rotor, is presented to provide the 
base for discussing the obtained results. Also presented in Chapter II, is a quick survey of 
the most relevart rererences who developed and used panel methods to study unsteady 


aerodynamics of oscillating airfoils in an incompressible, inviscid flow. 
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In the third chapter. the propulsive force is estimated according to the OH-6A helicop- 
ter parameters obtained from reference [36]. In this analysis, the blade is divided into sec- 
tions. Assuming hovering conditions. the reduced frequency at each section was evaluated. 
Using the unsteady panel code the time history of the airfoil drag ( or propulsive force ) 
coefficient was recorded. For the assumed plunging motion with plunging amplitude as 
obtained from reference {1], for a given mode of vibration, the average propulsive force 
for each section was calculated. 

The total average propulsive force was obtained by adding together the contnbution 
from each blade section. The obtained results showed that the resulting propulsive force 
from a single airfoil oscillating in plunge is not enough to explain the significant power 
reduction on the OH-6A helicopter. This led us to study the aerodynamic interference 
between two blades, which is the main subject of the fourth chapter. 

In Chapter IV, the two-airfoil unsteady panel code developed by Pang [ 25 ], 1s used 
to study the wake interference between two airfoils. This code has been formulated to 
solve for the potential flow for two airfoils executing unsteady motion in an inviscid 
incompressible flow medium. It is an extension of the code, developed by Teng [24], for 
single airfoils. The technique uses the well known panel method for steady flow and 
extends it to unsteady flow by accounting for the continuous shedding of vortices into the 


trailing wake. 
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Numerous case-runs are presented to illustrate the aerodynamic interterence between two 
oscillating airfoils. Results obtained for plunging and pitching airfoils show a significant | 
production of propulsive force at certain wake phase angles. In Chapter V the question of 
propulsive efficiency due to plunging or pitching and the possibility of single-degree-of- 
freedom pitching flutter are discussed. Chapter VI discusses the findings of this research 


and presents the main conclusions. 
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II. BACKGROUND 


The etfect of plunge and pitch oscillations on the resultant drag of an airfoil 
has been the object of aeroelastic research back almost to the time of the Wright 
brothers (1903). In 1909 and 1912, respectivelv, Knoller [10], and Betz [11] were 
the first ones to explain the birds’ ability to generate a forward thrust by means 
of wing flapping! Another one of the earliest investigators of this era was R. 
Katzmavr [7], who in 1922 conducted two senes of experiments to explore the 
etfects of oscillatory fluid or airfoil motion on resultant lift and drag. 

Katzmayr’s experiments which were performed in the aerodynamic laboratory 
of the Technical University of Vienna, may be divided into two series. In the first 
series, the angle of attack of the wing model was changed by oscillating the wing 
model about an axis parallel to the span and at right angles to the airflow. In the 
second series of experiments, the airfoils were stationary while the airflow itself 
was subjected to periodic oscillations. Here. for certain air velocities and frequen- 
cies, and considering thick airfoil sections, Katzmayr measured negative drag or 
positive propulsive force. This discovery, to be later known as "Katzmayr effect". 


helped to explain the propulsive mechanism of birds flying through the air. 


l Platzer. M. F.. Neace.K. S.. and Pang, C. K.,"Aerodvnamie Analvsis of Flapping Wing 
Propulsion", AIAA paper No. 93-0484. 31st Aerospace Sciences meeting, Reno. NV, January 


11-14, 1993. 
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Later. E. G. Richardson [8] used the "Katzmavr etfect" to explain thelocomotion 
of fish through water. In his paper, Katzmayr summarized his results to show: (1) 
the effect of flowing air, whose direction is undergoing constant periodical change, 
iS quite favorable. (2) the wing section which exhibits favorable characteristics in a 
constant airflow. works better in an oscillating flow. 

In 1935, Theodore Theodorsen [3] introduced the general theory ot aerodv- 
namic instability and the mechanism of flutter. In the first part of his paper. 
Theodorsen developed the velocity potentials due to flow around oscillating air- 
foil-ailerons. In the later part of the paper, he developed the differential equations 
of motions and determined the flutter speed. 

The most relevant work with respect to the present research is that by I. E. 
Garrick [5], who applied the work of Theodorsen [3] to obtain a closed form ana- 
lytical expression for the propulsive force generated by an oscillating airfoil. An 
airfoil oscillating in pure plunge or flapping motion was shown by Garrick to have 
a propulsive force throughout the entire range of reduced, k, while for an airfoil 
oscillating in pitch, propulsive force can only be achieved for a specified range of 


reduced frequency, k, where k must be greater than 2.0. 
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In 1942, W. Schmidt [14] recognized that the Knoller-Betz or Katzmavr etfect 
applies to both a flapping airfoil in uniform flow or a stationary airfoil in oscillat- 
ing flow. He demonstrated experimentally that a stationary airfoil positioned in the 
wake of a flapping airfoil increases the propulsive efficiency to almost 100 percent 
throughout the whole frequency range because the stationary airfoil converts the 
vortical energy venerated by the flapping airfoil into additional thrust. He also 
found that the mechanical limitations inherent in pure flapping motions could be 
overcome by an arrangement which he called the wave propeller. 

H. Bosch [15] investigated the interference of two lifting surfaces in two- 
dimensional unsteady incompressible flow, using classical flat -plate theory. He 
showed that a harmonically flapping airfoil upstream of a stationary airfoil increases 
the propulsive efficiency to almost 100 %, thus confirming Schmidt’s expenmental 
findings. Bosch’s analysis is limited to flat-plate airfoils oscillating about a zero 
incidence mean position. 

In 1957, Loewy [4] developed a two-dimensional model for the unsteady 
aerodynamics of the blade of a hovering rotor. He derived an expression for the 
two-dimensional oscillatory lift deficiency function as a function of the reduced 
frequency, wake spacing and the ratio of the airfoil frequency of oscillation and the 


rotational speed. 
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The wake in hover or vertical flight consists of helical vortex sheets below the 
disk, one from each blade. Unsteadv motion of the rotor blade will produce shed 
vorticity in the wake spirals. | 

Figures (2.1A) through Figure (2.1D), from Ref. [4]. Loewv illustrated the‏ کا 
basic elements of the unsteady flow field of a helicopter rotor. A rotor in hover-‏ 
ing or vertical flight produces a trailing tip vortex with a downward axial veloc-‏ 
ity that, if otherwise undisturbed, forms a contracting helix as shown in Figure‏ 
(2.1A). If the inflow over the rotor disk is constant, then the fluid off the trailing‏ 
edge of the blades makes a helical surface with horizontal radial elements as shown‏ 
in Figure (2.1B).‏ 

If there is an oscillation in blade effective angle of attack, blade lift will alter- 
nate also, and as a result of these changes in lift, vortices will be shed continu- 
ously at the blade trailing edge. These vortices fall along the horizontal radial 
elements of the helical surface shown in Figure (2.1B), so long as the oscillations 
in angle of attack are small. Figure (2.1C) illustrates this local sheet of shed vor- 
ticity. It should be noted that vorticity is considered to be on the helical surface 
shown in Figure (2.1B); the vertical displacements from that surface shown in Fig- 
ure (2.1C) represent the strength of the vorticitv at a particular azimuthal and radial 


position [4]. 
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Figure (2.1) Loewy's Schematic. Elements of Unsteady Rotor Flow Field 


1.2. 


The variation in this vertical displacement (hence vortex strength) around the 
azimuth corresponds to the historv of the motion of a given blade element at a 
fixed radius. Variation of shed vortex strength in the radial direction at any fixed 
azimuth angle is a function of the variation with blade span of (1) blade planform; 
(2) amplitude of oscillation of effective angle of attack, and (3) relative air veloc- 
ity. Since (from the Helmholtz theorem ) vorticitv cannot begin or end in space, a 
variation of shed vorticity in the radial direction implies that there are trailing vor- 
tices at constant radii similar to and inboard of the tip vortex. These trailing vorti- 


ces have been induced in Figure (2.1D). 


20 


The schemauc drawings, Figures (2.1A) through Figure (2.1D), Ret. [4], indicate 
the difficultv of attempting to obtain a complete representation of unsteady rotor 
aerodynamics. One means of simplifying the problem is to consider extreme "pitch" 
values for the helix. For the case of very high inflow velocity, u, in relation to Q, 


and the opposite, where u is verv low compared to © This is shown in Figure 


When the vertical spacing between adjacent helicai surtaces ot shed vorticity is 
very large, then we expect that all shed vorticity beyond a small fraction of a rev- 
olution would be too far below the blade in question to have a significant effect. 
Under these conditions, 1t would be sufficient to account for only the attached 
vortex sheet within that fraction of a revolution, as in Figure (2.2A). 

On the other hand, when the vertical spacing between the adjacent helical 
surfaces of shed vorticity is very small, all the sheets of shed vorticity tend to pile 
up on each other. and the effect of that vorticity close to the blade in question 
( shed by blade passes ) 1s of more importance than that which exists beyond 
a reference azimuth angle on either side of the blade. This situation 1s depicted in 


Figure (2.2B), Ref. [4]. 
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A. High Inflow B. Low Inflow 


ligure (2.2) Schematic Representation of Unsteady Rotor Flow Field. 
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With. low rotor thrust values ( low collective pitch ) the wake remains near 
the rotor disk and therefore passes close to the following blades. Thus. the shed 
vorticity sheets below the rotor disk must be accounted for to correctly estimate the 
unsteadv loads. For high inflow or high-speed forward flight the rotor wake will 
be convected away from the blades. The wake far trom the blade section will have 
little influence. so emphasis is on modeling the wake near the blade. In the case 
of low inflow, this consists of vortex sheets that are nearly pianar surfaces paral- 
lel to the disk plane [4]. 

Based on these considerations. a two-dimensional model for the unsteady aero- 
dynamics of the rotor can be constructed. It is assumed that the chord, amplitude 
of the oscillation in effective angle of attack, and relative airspeed vary slowly 
enough with span so that what occurs aerodynamically at one blade radius station 
is essentially duplicated on either side of it. 

The assumption means that the flow problem at a given blade section is two- 
dimensional. Furthermore, consistent with the idea that only the vorticity near the 
blade section has an important effect, one may allow the planar rows of vorticity to 
extend to infinity in the horizontal direction in order to achieve mathematical sim- 
plification. It should be stated that this is a first-order theory, so that such effects 
as that of the wake upon itself have been neglected. The oscillatory flow induced 
at the airfoil in the free-steam direction, for example, will affect the airfoil lift and 


moment as well as the vorticity shed at the trailing edge. 
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Since the shed vorticity is being accounted for as perturbations in the main 
flow. this component of induced velocity will be small compared to the free-steam 
velocity (except when the spacing approaches zero) and must react with small. first- 
order airfoil displacements. Because of the physical unrealities associated with the 
standard mathematical model representation of the problem (such as zero thickness 
of wake and airfoil, zero displacement of the wake from horizontal planes, and 
infinite velocities at the core of a vortex), the accuracy of results as the wake 
spacing goes to zero is questionable. The effect of the horizontal component of 
induced flow. therefore, has been neglected in. Loewv's analysis. 

Shown in Figure (2.3) is the two-dimensionalized model of unsteadv rotor flow 
[4]. The portion of the circular surface which is determined by (1) a particular 
blade radius, (2) the azimuth angle on either side of the blade section (within 
which the shed vorticity is of importance), and (3) the vertical T spanned by 
a given numbers of rows of vorticity, can be regarded as a plane, one in which 
the two-dimensional unsteady aerodynamic problem may be attacked. 

This is shown in Figure (2.3 A), with. an arrow indicating that the rows of 
vorticity under the rotor disc extend to infinity. Consistent with the idea that only 
the vorticity near the blade section has an important effect, one may allow the pla- 
narrows of vorticity to extend to infinity in the horizontal direction in order to 
achieve mathematical simplification. These final modifications are shown in Figure 
(2.3 B), which is a two-dimensional model, Ref. [4], of unsteady rotor aerodynamics 


for a single-bladed rotor operating at low inflow. 
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Figure (2.3) Two-dimensional Models of Unsteady Rotor Flow. 
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The general approach used bv Loewv is the same as that of Schwarz [30]. 
Only the etfects resulting from oscillatory motion of a thin airfoil are considered, 
and the airfoil itself is represented by continuous vorticity distributed along a 
straight line [9]. It may be also noted that the effect of viscosity in damping the 
shed vorticity could be accounted for by considenng only a finite number of rows 
beneath the reterence airfoil, or by multiplying the vortex strength by an experi- 
mental decay function related to the time through the rotational speed and the 


number of rotational cycles corresponding to a particular vortex row. 
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This has been omitted from the analysis since vortex decay time is long 
compared to the rotational periods of most rotors. An extension to Loewy 's work is 
that by Shipman and Wood [29] in 1970. The authors present a method for 
determining rotor blade flutter in torward flight. Incompressible unsteady aerody- 
namic theory was applied where shed vorticity both in and below the plane of 
the rotor is accounted for. This was made possible by noting that the damping for 
rotor blade aeroelastic instability vanes with velocity. Since the tangential velocity 
of any blade segment varies with azimuth. so will flutter damping’. 

It 1s assumed at the onset of flutter that oscillations will begin to build up 
prior to the blade reaching a critical azimuth position, then decay as the blade 
moves beyond this point. This build up and decay means that the vorticity shed 
due to the oscillations will be contained within a double azimuth region on either 
side of the critical azimuth position. Assuming this region to be small allows the 
wake system to be two-dimensional. 

The lift deficiency function resulting from the two-dimensional wake model is 
compared with earlier results obtained for a helicopter in hover and fixed wings. In 
the limiting case when the advance ratio is very small, Shipman and Wood's lift 
deficiency reduces to Loewy's lift deficiency function. Also in the limit as the 


wake spacing is very large, this function reduces to Theodorsen's lift deficiency 


function . 


l Shipman, K. W. , and E. R. Wood" A Two-Dimensional Theory for Rotor Blade Flutter 
in Forward Flight." Journal of Aircratt, Vol. No. 8, Number 12. pp. 1008-1015, 


December 1971. 
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The theory was applied to bending torsion flutter for the tip segment of a rotor 
blade. The results followed the normal trends of having the flutter velocity decrease 
as the blade center of gravity was moved aft and of having the flutter velocity 
increase as the blade stiffness increased . 

The previously shed wakes are destabilizing so that they reduce the flutter 
velocity. The buildup and decay produces wakes that are essentially centers of vor- 
ticity. Thus their position with respect the reference blade will be sensitive to the 
advance ratio which determines the vertical spacing. It was found that as the 
wakes were brought closer by reducing the advance ratio and/or intlow ratio, the 
flutter speed is decreased [29]. 

The phenomenon of airfoil flutter is caused bv energy extraction out of the 
wind stream by means of airfoils which are freeto oscillate in pitch. and. plunge. 
As described by Duncan [19], this effect can be easily demonstrated in a wind tun- 
nel if a device is built which allows to vary the phase relationship between pitch 
and plunge. Phase angles near 90 degrees produce maximum energy extraction so 
that the device operates as aflutter engine. McKinnev and DeLaurier [20] demon- 
strated that this effect can be used to convert wind energy into mechanical energy. 
Scherer, McKinney and Delaurier used Theodorsen's flat-plate theory for perfor- 
mance estimates of the flapping foil propulsor and of the wingmill, respectively. 

It was shown by W. Schmidt, that optimum energy extraction occurs for har- 
monically flapping airfoils upstream of a stationary airfoil. There is a need for 
aerodynamic analysis methods which enable systematic studies of the effect of air- 
foil geometry, incidence angle variation, amplitude of oscillation and airfoil interac- 


tion on the aerodynamic forces due to incompressible flow over airfoils or airfoil 
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combinations executing pitch or plunge motions. A panel method was developed in 
order to do parametric studies of these effects. In the past, a number of investiga- 
tors have solved the steady flow problem using source and vortex paneling, the 
most prominent ones being Hess and Smith [21]. A few authors have extended this 
approach to the case of unsteady motion of single airfoils. notably Basu and Han- 
cock [22] and Kim and Mook [23]. 

Basu and Hancock [22] presented a numerical method for the calculation of 
the pressure distribution. force. and moments on a two-dimensional airfoil undergo- 
ing an arbitrary unsteady motion in an inviscid incompressible flow. The method 
was applied to (1) a sudden change in airfoil incidence, (11) an airfoil! oscillating in 
high frequency and (iii) an airfoil passing through a sharp-edged gust. 

Teng [24], developed a computer code for the numerical solution of unsteady 
inviscid incompressible flow over an airfoil. Pang [25] extended the work of Teng 
and developed a computer code to study unsteady airfoil interference effects. The 
technique uses the well known panel methods for steady flow and extends it to 
unsteady flow by accounting for the continuous shedding of vortices into the trailing 


wake. 
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Ill. DRAG REDUCTION DUE TO 


BLADE FLAPPING 


(KATZMAYR EFFECT ) 
A. INTRODUCTION 


In this chapter. a discussion will first be presented of the physics ot flapping 
airtoils which explains the bird’s ability to generate a forward thrust by means of 
wing flapping. This will be followed by a description of two methods capable of 
predicting the pressure distributions, lift, pitching moment, and drag or thrust on airfoils 


which execute harmonic oscillations in plunge or pitch. 


The first method is based on the generalization of the well-known steady-state 
panel method to the case of unsteady airfoil motion. It is based on the assumption 
of inviscid, incompressible flow. The second method is based on the classical The- 
odorsen theory for oscillating flat-plate airfoils in incompressible inviscid flow. The 
unsteady panel method permits the investigation of finite-thickness airfoils of arbitrary 
geometry which may execute rather general unsteady motions (ramp and oscillatory ). 
In contrast, Theodorsen's analysis is restricted to thin airfoils ( flat-plates ) and to small 
amplitude oscillations. The chapter ends with a comparison of results computed by both 


methods. 
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B. THE PHYSICS OF FLAPPING AIRFOILS IN LOW-SPEED FLIGIIT 


The phenomenon of thrust generation due to wing flapping or due to steady flight 
in a sinusoidal gust can be explained in a rather elementary way. Consider the latter 
case of flight through a sinusoidal gust. as shown in Figure (3.1). The relative wind 
vector changes direction between a maximum positive and negative angle of attack. The 
resultant aerodynamic lift is perpendicular to the instantaneous wind vector. If we as- 
sume the viscous drag to be small. then it is easily seen that the sinusoidal gust gen- 


erates a sinusoidally varying thrust. 


A similar effect is generated if the airfoil executes a sinusoidally varying plunge 
motion ( wing flapping ) about an otherwise steady flight conditions. This explanation 
of thrust generation due to wing flapping or due to flight through a sinusoidal gust was 
first advanced by Knoller [10] and Betz [11]. Its experimental demonstration was 
accomplished by Katzmayr [7] and is theretore generally referred to as Katzmayr 


effect £. 


Since it 15 well known that every change in airfoil incidence causes the shedding 
of a starting vortex, a more precise explanation of the Katzmayr effect has to take into 
account the continuous vortex shedding due to wing flapping or due to flight through 


a sinusoidal gust. 


l Platzer. M. F., Neace. K. S., and Pang, C.K.,"Aerodvnamic Analvsis of Flapping Wing Propulsion’. 


page 1, AIAA paper No. 93-0484, 31st Aerospace Sciences meeting. Reno. NV. January 11-14, 1993. 
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Figure (3.1) Propulsive force on plunging airfoil 


32 


Similarly, the quantitative prediction of the induced oscillatory forces (lift and 
thrust) had to await the development of a complete unsteady aerodynamic theory. 
Such a theory was first accomplished in the 1930’ s by Theodorsen [3] in the United 
States and by Kussner {16] in Germany. It can be found in the textbooks on aeroelas- 
ticity, for example the book bv Bisplinghoff. Ashley, Halfman [41]. Therefore only 
the final formulas will be given here and the reader is referred to those texts for 
additional details. More recently, a second approach to the determination of the 
unsteady aerodvnamic forces on oscillating airfoils became feasible with the advent of 
high speed computers and the development of efficient numericai solution methods. The 


fundamental building blocks of this method are presented in the next two sections. 
C. PANEL CODE FOR STEADY INCOMPRESSIBLE INVISCID FLOW 


In potential flow theory, the flow field around an airfoil may be represented by 
the velocity potential. Considering contributions from the free stream flow and the 


source and vorticity distribution, the total potential may be constructed: 
D = Op + Oz + dy (3.1) 


where 


US costas ss vs aA) e] 


Dp 





p = pe In (r) ds (3.2) 
TT 


5 ^ 


pass -[ O gas 
: on 
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The source distributions q(s) vary from panel to panel. while the vorticity strength 
7(s) 1s assumed constant for all panels. The value of representing the flow past an 
airfoil by surface singularity distributions lies in the fact that these singularity distri- 
butions automaticaily satisfv Laplace's equation. the governing flow equation for invis- 


cid incompressible tlow: 
DEO 0 (3.20 


Since Laplace's equation is a linear homogeneous second order partial differential 
equation. the superposition principle used in Equation (3.1) holds. The boundary con- 
ditions include flow tangency at control points (mid-points of panels) and the Kutta 
condition at the trailing edge, requiring equal tangential velocities for the first and last 
panels. By evaluating the integrals along the airfoil surface, the potential may be de- 
termined at any point in the flow field. Each point is defined at a radius (r) and angle 


(0) from a chosen reference point on the airfoil. 


The reference point in this studv 1s the leading edge. The airfoil is represented 
by a number of defined points. called nodes [34]. More points produce greater resolution 
and better accuracy. One hundred to two hundred points are usually sufficient, with 
the larger number used for more complicated airfoil shapes or more involved calcula- 
tions. The lines connecting these nodes are the panels. There are ntl nodes with first 


and last node overlapping. 


34 


"gure (3.2) depicts the panel geometry. Numbering starts at the tratling edge. then 
resses along the lower surface. leading edge, and upper surface, and ends at the 
ng edge. The unit normal vectors (n;) are perpendicular to the panels and directed 
ard from the atrforl surface. The unit tangential vectors (t,) are parallel to the 

Is and the positive direction is defined with increasing numbering {n to ntl) 
panels may vary in length, with the exception of the first and last panels, which 


be equal im order to use the Kutta condition at the trailing edge. 





Figure 3.2 Panel Method Geometry 
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The required input consists of the number of nodes on the aırloli surface. the 
coordinates or the nodes referenced trom the ieading edge. and the angle or attack in 
degrees. The program produces normalized velocities and pressure coerficients at each 
control point as output. The use of intluence coeitficients leads to a. straightforward 


procedure tor programming the equations. 


An aerodvnamic iniluence coetlicient is. detined as the \velocitv induced at a fieid 
point by a unit strength singularity distnbution on one panel. Fer the two-dimensional 


steady ¿low problem. the tollowing tntluence coefficients are needed [341 


A". normal velocity component induced at the i" panel control point bv unit 


strength source distribution on the ie panel 


| | : DE "D 
A” = sin (0. - 89.) In + cos (9:- 09 تت‎ . G (3.4) 
ا ا‎ V P ٢و‎ ٢۳٢ 
l 1] | 
= 1/2 | != | 
اد‎ : tangential velocity component induced at the n panel control point bv unit 
strength source distribution on the p panel 
[ 1 f m + | 2 5 
A.. = — sin(0.-0)p..—-cos(8.-0.))]In—— , 14) (5:59 
t a 1 ] 1] 1 | 5 
A i: Xs 


. . . . . 
B”: normal velocity component induced at the i" panel control point by unit 


strength vorticity distnbution on the T panel 


| A 
= 2 cos (6) 20) in = sin Gt EE جم 1 :كح‎ (3.6) 
1 1 ] 1 ] 1] 


I 
)ا‎ 
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: | ፡ 
Br tangential velocity component induced at the 1" panel control point by unit 


strength vorticity distmbution on the Da panel 


|| 


| | EI | 
B. cos (O — 0 ) D.. - sin(0O. —9 )In d 
1] - | ይፈ |] | | 


=f )( 
2 e 
B ا‎ 
2چ‎ m 
where the veometncal quantuties. depicted in. Figure (3.3) are detined bv 
B = Na کا‎ 
2 | || 
c= (xm, - y Pans = ኪረ: e 
N m \ S | کا‎ 
ےہ ہت‎ rd c 20 
vm. = —- 0. = arctan tw) 
: 8 | N -x 
ከ 1 
/ t f =y. 
| pt y; * | 2 لا‎ 
B See tall) == | arctan | — 
MX J wm; - x, 
The first boundarv condition requires flow tangencv at the control points 
1 O ET n (3.9) 
In terms of intluence coetficients (with Pe = lo) 
n n 
n 1 : : 
A..q. +v N° B.. - -0.) = =1,2. n x | 3.10 
> | iji. DI Bi, > sin (a 0.) TEE 129 ۱ (3.10) 
ا ا‎ 


The second boundary condition is the Kutta condition. which states that the pressures 


on the lower and upper panels at the trailing edge must be equal if the flow is to leave 


the trailing edge smoothly. 
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Figure (3.3) Relationship Between Geometrical Quantities 
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Using. Bernoulli's equation; the pressure coetficient can be expressed as 


D 
D a D 
۸ | eMal 
uH 03 
P | 2 31 F. 
DV Y^ 


als) 


The pressure equilibrium also implies equal velocities for incompressible flow. Since 


the noimal velocities are taken to be zero, the boundary condition may now be stated as | 
rl A [ 7 
(V y = Ne) (20129 


where the negative sign ts strictly due to the adopted convention of positive tangential 
velocities in the direction. of increasing node numbering. Since the flow is positive to 
the tight, the panel downstream of the front stagnation point will have negative values 
for computational purposes only. It is important to note that not all the lower surface 
panels have a reversed sign, only those downstream from the stagnation point. This is 
especially significant for non- symmetrical airfoils or any airfoil at an angle of attack. 


In terms of influence coefficients, the normalized equation becomes 


n n n n 
6፡1 ip ا د‎ 
- E | 1 | ፦ኘ > Bi; - cos (a. — 0) = E | د غ‎ E > B, 605 (a - 0,) (3 13) 


0۴ E = rek 


Equations 3.10 and 3.13 represent a linear algebraic system of (n+l) equation and 
(nt 1) unknowns. The unknowns are the source strengths which vary from panel to 


panel ( q qQ4) and the vorticity strength y. 
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Expanding and rearranging Equauon 5.10 for an exampie airroti or n —772 nodes and 


- 


saneis results in 


ር 11 1 نے‎ Am - 

Avus c rece MEE SB ር ርን የከ حل‎ EJ ከከ. 

1 4 lee yan =! IM = an 

a in. | 1 p ወው. ሰበ 3 
E As de لو‎ በመከ ፕ ከርሮ TD c = IES ^23) [ጋ 14) 
2 il il 4 ad 
AI. ووا‎ ٤ ټم‎ o UB Be ۱۱۱ ም ሞ ሺር ፡ች፡፤] EIE ይ 


The equauons now readily lend themselves to solution in matrix torm Recasting 
‘sith simpler notation. the A", terms (coefficients or a, 1 mav be renamed a., and the 
sum of all. B" ,, terms in parentheses (coefficients of + J renamed a,,.;. where i = |. 
2 ..n and !21.2. . ..n. The terms on the right side ot the equation may 


be renamed b,. The (nti) equation, or in. this example. the pu equation comes 


from Equation. 5.13 in a similar manner: 


( 


t l i ‘ 
8] ۰ 1) تب نلک‎ 0 


p = E 


1 ( 0 5'-ን ; ከ ( 5 |13ፐ Dee TEE E ( BT B! =, ~, ) | 
= co (a -O ة7 ہہ‎ (3.15) 


The coetficient q; mav be renamed a-,, All ot the B' terms in the brackets together 
torm the coefficient of y, now renamed a-,7,. The entire right side of the equation 


constitutes the new term b-,. 
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Finallv expressing this system can be written Tn concise. matrix. form 


E es Sd] ٩ 5i 
NIME ae R Da 
Lx E 000 ET (3-16) 
EM LS s 9, b, 
EIA E 0ط‎ 7 Dess] 


This system is solved in the program using a Gaussian Elimination subroutine. With 
the values ot the q, and y known, the velocity at each panel control point may be 


calculated 


n 
ከ l ee : => 1 = 3 0 
V. 29) a A A qme (3.17) 
دا‎ 
The total velocity is equal to the tangential velocity due to taking the normal velocity 


to be zero. 
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D. UNSTEADY PANEL CODE 


The extension of this method t5 an urtoil expenencine 9 rae ካጣ) 
ingle o; attack can be accompiisned as toilows. it 1s important to recall Pranatl’s 
eariv tlow ‘\isuaitzation experiments “nich demonstrated the senerauopng በከ 
starting vortex orf the irailing-edge whenever the angle Cr attack was changed. 
Hence continuous changes in angle or attack produce 2 conunuous snedding of 
vorticity into the trailing wake. This shedding can be exptainea by tne Heimholtz 
vortex theorem which requires that anv change tn the circulation around tne airtoll 
must be matched bv the appearance ot an equal counter vortex ar the trailing 


edge in order to achieve constancy or the total circulation in tne iON 


The vortices shed from the trailing edge move with the tluid particles of the 
surrounding tluid and hence are swept downstream with a speed essenuallv equal 
to the free-stream speed. Therefore these vortices will stav close enough to the 
airfoil for a finite time to influence its pressure distribution. Any unsteadv airtoil 
theory theretore must describe this vortex shedding process. It 1s this feature 
whicn distinguishes unsteady airtoil theory trom its steady counterpart’. The 
above considerations suggest that the airfoil can be modeled bv similar source 


and vortex distributions as in the steady case. 


Platzer. M. F. Class Lecture Notes. Naval Postgraduate School. Monterev. 


CMOS ODD 
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Hence source panels are again placed on the airtoii surtace together with a 
single vortex distribution which is the same for each pane: ¡An arbitrarv change in 
angle of attack then is subdivided into small step changes such that a starting 
vortex is shed into the wake each time a step change occurs and the tlow tangencv 
and Kutta conditions are entorced at each time step. 

In addition to the n airtoil panels a wake panel 1s assumed to be attached to 
the trailing edge such that. atter each time step. the vorticity of the wake panel is 
assumed to be concentrated into a single point vortex which detaches trom the 
airfoil and starts moving downstream with the free-stream velocity. As one wake 
panel detaches from the airfoil a new one is created and a wake of point vortices 
1s created. As in the steady problem n unknown source strengths and one unknown 
vorticity strength on the airfoil are introduced. 

In addition, there is now the unknown vorticity strength of the wake panel 
whose length and onentation are also unknown. Hence there are three additional 
unknowns which need to be determined by a proper set of equations: (a) the 
flow tangency conditions provides n equations. (b) the Kutta condition provides an 
additional equation, (c) the Helmholtz theorem provides an equation tor the vortic- 
ity strength on the wake panel. 

Therefore two more equations are needed to determine the length and orienta- 
tion of the wake panel. These are obtained by making the following assumptions. 1.e.. 


(d) the wake panel is oriented in the direction of the local resultant velocity at its 
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፥ 


nid-pomt. (e) the length of the wake panei is proportrona! to ihe magnitude ul the 


.ocal resultant velocity at its mid-point and to the ume-step. The Kutta conaition 


and the last two conditions are nonlinear equations. and theretore necessitate an 


i22]. Shows the panel method 


| - 


iteratrve soiutton procedure. Figure (3.4). from Ret. | 


representation for unsteadv row. 


44 


Vortex Shedding at [ime Step t, 


Helmholtz s theorem 
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Figure (3.4) Panel Methods Representation for Unsteady Flow 
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٠٠.) SIE 
As stated berore. Theoaorsen s theory tor oscitlating jiat-ptate airoils in inviscid 


ncompressibcle tow is well documented. ‘hererore. | e summanze rere o niv [ከይ 


tinal tormuias. The ntt on an airfoil oscillating at treauency © ts given Cv 


1 


፡ክ)ህ፡.፡.ን፡.፡. 1! 11፤.. ፡ | (2.48)‏ اط نونج سنا 
ec‏ 

ሸህ. ም 

ان اکاک کا ل Li. = DI Go‏ 

k= ob/U «q bis the semi-chord and "a" 15 the elastic axicm 

C(k) is the lift deficiency function and is given by 

Hi^ (k) 
C(kK) 9 ————— ——— = F(k) - 6 (ا)‎ (3.19) 


D 


1.1) 


where H“ tk) is the Hankel function of second kind. and is given in. terms 


of Bessei functions as. 080 =) JS) 9 EUN 
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F. GARRICK'S PROPULSIVE FORCE ANALYSIS 


Garrick [5] applied Theodorsen s theorv to determine ihe thrust and drag 
generated bv an airioil which executes a sinusoidal motion in pitch or in. piunge. 


He found the tollowing formula for the average propulsive force 
P= robo” 7ہ‎ ~ (a, +b.) ae —? (a, Di) Calle 15290) 
where e Eos 
NE ۶" ١٣ یپ۶۹‎ 7 
MED GIL CNC O UC ACOSO OF 
uu casam 9  - + 70 
BE Dal BE 200 ۱مم رہ-‎ [ 
by = (b/2)[ (I2 + G/k ) cos ( زىم حدم‎ ^ (F/K) sin ( p5 — رہم‎ 
For pure plunge, the average propulsive force is reduced to 
P. - xpbo^h. F'(k) «G (k). (3.21) 
For pure pitch, the expression for the average propulsive torce is given bv : 


ን ም 
D up» D cua. (3022) 


where a; and b» are defined in Equation (3.20). 
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ihe expression tor the average drag in Equation (35.20) «t Garnck s report 
3:07 [5]was tound to be : error The correct expression as given ከ 


O ጀ O AS 


0۰٦ | = i ] r | Gi 
MES C A EN = -F Eg t «s TR i M D) 


و 
S‏ 
رو 
» 
J‏ 
= 


) E =>... Tc DNE, 
G. COMPARISON OF RESULTS 


ín a recent thesis. Riester [ 43 ] presented a detailed comparison ‘vith 
Theodorsen s predictions ot the lift and pitching moment computed with the un- 
steadv panel code. Generally. good agreement was found. thus lending credence 
to the reliability of the panel code. In the present work. we are interested in 
the prediction of the thrust and drag forces generated by an oscillating airfoil. 
Calculations were performed to ascertain the panel code’s capability. Figures 
(3.5) to (3.7) show the companson between Garrick s predictions and the vanel 
code results. These figures present plots of the average propulsive torce coef- 
ficient vs the non-dimensional plunging amplitude of oscillation. HBAR (HBAR 
= h/c ), for three reduced frequencies, k=0.1. ! and 2. The two methods are 
in good agreement for small plunge amplitudes. As might be expected. the 
agreement starts to deteriorate with increasing plunge amplitude. ¡he panei code 


was applied to an NACA 0009 airfoil using 100 panels. 
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Figure (3.5) Comparison of Panel Code Results with Garrick ( k=0.1 ) 
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Figure (3.6) Comparison of Panel Code Results with Garrick ( k= | ) 
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Figure (3.7) Comparison of Panel Code Results with Garrick ( k=2 ) 
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H. ESTIMATED POWER REDUCTION DUE TO HHC FOR OH-6A 


HELICOPTER 


The ultimate objective of this research was to determine the erfect or unsteady 
aerodynamics on the measured power values for tne OH-6A helicopter. “Ve are 
especially interested in knowing whether these etrects could produce the power 
changes recorded dunng the NASA - Army open loop ilight tests in the. light 


regime trom hover to 100 knots. 


In order to estimate the change in torque due to propulsive force ( zoot- 
pounds per blade) due to HHC. the blade span was divided into sections As- 
suming hover conditions. the reduced frequency at each blade section was cal- 
culated The panel code was used to compute the propulsive force coefficient. 
Cd, at each blade section assuming pure plunge motion with a plunging am- 
plitude of 1". The blade deflections. obtained from Wood et al [27], were then 


used to estimate the real value of the plunging amplitudes at each blade segment. 


Finally, the net blade torque could be calculated by considering the propulsive 
iorce at each blade section and multiplying this value umes the rotor radius. then 
summing over the length of the blade. Table 3.1. from Ref. [42], shows the 
OH-6A helicopter parameters that were used to estimate the propulsive forces. 
The most important parameters are: number of blades. Q =4. rotor radius (ft). 
R= 13.2. blade chord (ft), c= 2b = 0.57. rotor speed (rad/sec), £2 لل لل‎ 


helicopter weight (lb), G= 2550. 
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TABLE 3.1 THE OH-6A HELICOPTER ROTOR PARAMETERS : 


رما 
)1 


Rotor radius (ft) | 


Number ot blades 4 
Rotor speed (rad/sec) 3096 
Blade chord (tt) ه7‎ 
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Table (3 2) shows, at each blade section r ከዚ. the reduced ٤+ الل‎ 


local speed. Ur), the computed average propulsive rerce. P..ir) ana the torque. 


+ م 


inu) Jue to tha propulsive torce at each 3۱806 <1: - - ٢ 
nated torque at each blade section. :he resuitant torque due to : ۹ ٥ ٦ 


O 4 


oer blade ‘was round :0 ከ6 8 11 ft-lb assuming « pure viuneln: 


ہے 


mouon with 


(y 


۱" 


constant plunging amplitude or | 


.-- 
مد 


"teach blade )۰ء۰‎ +.٦ 
plied Garrick’s equation. and considered the etfect or the second to tweltth har- 
monic loads. He tound. ror 2P harmonic ( HHC " on ^ —20 decree Sams ME 
Knots ). ihat the power gained is about !1.5 horsepower. Tor more cetari. the 


aa 


۲٥۵۹٢١٥٢٣٣۰٠٢ ٦7216۲160٣٢٢٠٠٠٢ + +++3+ؤ‎ 8 ۲ 


Ti: ^ O 4 5 2801 
1 ^ 


! ۵ا‎ to. 2 ! 1 now ٣ 
detlections tor several harmonics used by this reterence.. -nown on Figuret3.12) 
is the average drag coetficient. Cdav. versus reducea frequency as obtained with 
the panei code for a plunging amplitude ot 1". The airroil used was NACA 0015 


(OH-6A airfoil). The propulsive force. Pj, 1s given bv 


Bee Oel 7 (3.24) 


and the torque. Txh. at each blade section ts given bv: 
7 = Pos Ü E 


where U is the local speed at each section, C, is the propulsive torce coef- 
ficient as obtained by the unsteadv panel code and r 1S the radial distance of 
each blade section. For the OH-GA. the angular ireauencv . £2 = 5026 6 


therefore the power per blade is tound to be 0.72 horsepower as given bv 


PEZ A (3.26) 


o4 


TABLE 3.2 LOCAL REDUCED FREQUENCY. ik, (r). PROPULSIVE 


FORCE. Pq (r) AND TORQUE. Tq (r). 


r/R k (r/R) U(r/R) 


٢٢ (r/R) 








0 1 0 863 66 4 0.0618 0.122 
02 outs 1/5 132.8 0 0798 0 26 
0.3 0.2877 199 2 0.9530 0 44 
04 0.258 ^65.6 0 108 0 64 

| 05 TRIS EA 0.118 0.86 

| 06 0.144 398.4 0.127 1.09 

^ 07 0.1233 164.8 0. 15% 
0.8 0.1079 531.2 0.139 1.56 
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Figure (3.8) Blade Deflection vs Radial Station at 2P 
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Figure (3.9) Blade Deflection vs Radial Station at 3P 
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Figure (3.10) Blade Deflection vs Radial Station at 4P 
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FABLE 3.3 : Cd, versus REDUCED FREQUENCY د‎ at hbarz1" 
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IV. ANALYSIS OF WAKE 


¡NTERFERENCE EFFECTS 


A. INTRODUCTION 


The results ot Chapter HII show that propulsive torces. zenerated solelv 
by plunge oscillations of an airfoil are insufficient to explain te power 
reduction observed and measured on the OH-6A helicopter. Thererore. in this 
chapter. ve examine whether the power reduction might be caused 5v tavorable 


wake interference effects due to the wakes shed from preceding blades. 


To this end, the extension of the unsteady panel code to the case ot two 
airfoils [25] was adopted and applied to the study of the wake interference effects 
between two helicopter blades. Wake interference effects were studied in the 
1950"s by Loewy [4] by extending Theodorsen s analvsis. Loewy treats an 
infinite number ot layers ot shed vorticity placed beneath the rotor at a given 
wake spacing and developes a modified lift deficiency function. His interest 
was in the possibility of wake-induced blade flutter, which restricted his analvsis 
to the determination of lift and pitching moment changes due to wakes shed 
from preceding blades. A closer inspection of Loewy’s paper suggests that it 
should be possible to use the same approach to examine drag and propulsive 


force changes on the reterence blade due to wakes shed from neighbouring blades. 
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0111661121) cannas one Step further and investigate the 
case of a single wake shed trom the preceding blade which interteres with the 
wake shed from the reference blade! This case allows a direct comparison 


with the unsteadv two-toil panel code. 


Also. there are regimes of helicopter flight where oniv the tirst laver ot 
shed vorticitv is of consequence. In the following sections we tirst describe the 
two-foil panel code and Couch's modification of the toewv anaivsis. This 1s 


followed by a presentation and comparison ot the major results. 
B. UNSTEADY TWO-FOIL PANEL CODE 


The extension of the single airfoil code ( previousiv described ) to the case 
of two airfoils requires no new building blocks. However. the two-toil analvsis 
requires the introduction of five frames of reference, namely two moving local 
frames of reference which are attached to the two airfoils, two frozen local frames 
of reference. and the inertial frame of reference. Furthermore. it requires the 
sausfaction of the two Kutta conditions which are coupled non-linearly The 


solution requires an iterative procedure to compute the two vorticity distributions. 


I A detailed development of this theory 1s denved and presented 
in a recent NPS thesis bv Mark Couch. Ref. [40]. 
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It also necessitates the creation of a subroutine whicn transrorms ail coor- 
dinates in either or two local frames or reference to thé Global frame cr rererence. 
Finallv. it requires the eXtension of the influence coetticient concept to inciude 
ihe etfect or the second airfoil with its own wake and it requires the introduction 
ot an additional influence coefficient. namelv that on the wake element jue to 
the wake element from the other airfoil. If each airtoil .s modelled ah N 


panels. then this produces a 2Nx (€ *—5 ) matnx which ,s subsecuentb ) 


ov Gauss eliminauon 


Pang [25] was onlv able to venty this code ^v .omparne +٣٥ 
pressure distribuuons with. previous work bv. Giesinge 435] Hence tt was nec- 
essary to further evaluate the code bv companng its output with. other known 
solutions before applving it to the oscillatorv blade interterence problem. These 


comparisons are given in section C. 


The main parameters used in this code are: the airtoil tvpe. number of 
panels. N, the relative locations of the two airfoils. initial angle of attack. am- 
plitudes ot oscillation (pure pitch. pure plunge or both). irequencv ot oscillation 
for each airfoil, the pivot points (for pitch cases), rise timet for ramp motions ). 


These input parameters are in the file "fort. I" 


This program produces an extremely large amount of output to the screen. It 
IS convenient to write the output to a file on. a Unix. based. machine during 
program runs. The output file can get very large for long program runs (when 


the time historv is long). 
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For this reason. ihe logical vanable “output” was added to the input file in two 
modes ' true” or ' false". When this variable is set to ‘false’. the screen output 


is not printed and the required outputs are redirected to output tiles. 


The following list deseribes the output tiles and the data thev contain 
tort. 2. This is for user supplied airtoil coordinates. i: desired 
for. This tile contains the global coordinates ot AF1 at each time step. 
fort.+ This file contains the global coordinates of AF2 at each time step. 
tort.? This file contains the lift. moment, and drag coetficients tor poth airroils 
fort.S: This tile contains the pressure coetf. for AFI at each time step. 
tort.9° This tile contains the pressure coeff. for AF2 at each time step. 
tort.10: This tile contains the first airfoil s core vortex (wake) positions. 
fort.11: This file contains the second airfoil’s core vortex (wake) positions. 


fort.12: This file contains the angle of attack at each time step. 


fort.13. This tlle contains the y-axis translational motion at each time step. 


C. LOEWY'S ANALYSIS 


Loewv [4] developed a two dimensional model for the unsteadv aerodvnamics 
of the blades of a hovering rotor, including the effect of the shed wakes. Figure 
(41). trom reference [40], shows his two-dimensional model ot the helical 
wakes. The wake in the case of a multi-bladed rotor in hover consists of 


helical vortex sheets below the disk. one from each blade. 
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Assuming sinusoidal motion, the downwasn over ine rererence MO 


ITS vorticity (either bound or shed) mav be written. respectiveiv as 
NV = 3 S : “ያ E ع‎ | cf |) 


:where the barred quantities are complex. The induced velocity at a point x on 
the airfoil resulting from an element of vorticity of strength. 7, , 4t a. general 
point in the wake is given bv [4] 


Ay 
à 


m ۳ 
dw (<.t).™= 42) 


^ ^ ^ 


ከው ا‎ (70 ٣ 
where Q is the total number of blades and h is the verucal distance between 
successive rows of vorticitv. n is the number of revolution index. and q 1s the 
number of blades index as shown In Figure 4 |. The wake spacing, h 1s defined 
as [4]: 

zu 
J| bQQ 


where u is the intflow velocitv (time-averaged normal to the disk). 
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Figure (4.1) Aerodynamic Unsteady Model for a Multi-Bladed Rotor. 


The total downwash over the airfoil can be expressed conveniently by summing 
the integrals involving the bound vortcity and the row of vorticity in the plane 
of the reference blade with those for the rows of voruces below the plane of the 
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The first two terms on the right-hand side ot this integrai equauon are those 
that arise in the classical unsteady airtoil theory ort Theoagorsen | 3 |. the third 
and fourth terms contribute the influence or all the vorticity below the »iane 
of the rotor disc. Loewv solved the integral equation using Sohngen’s inversion 
formula [31]. and obtained his iift deficiencv function... C ( k.m.h). wnich ts given 


bv 


Ck inhi ee) a G a ae ee 2 o BER 


where K. is the reduced trequencv (detinea by 50/L. m is S15 NN 


rato Q., h is the wake spacing defined betore. and the real and imaeginarv 
parts of Loewy s lift deficiency function. F 9G are oven by [ie 


(Y, + BJ,)B‏ تكن ل 


P (E m ha . ። 
A -B 
(Y, - BJ) A - J,aB 
)وا‎ ٢ ٢٢ و‎ CET 1 
A -B 
where 
ا ا‎ eB Y APN 
and 
ا‎ E 
E ASIN (TIN) 
dE iL 7 . A . መጠ 
e" —2cos(2xzm) +e e COS Em cee 
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Here ıt can be seen that as the spacing between rows Of vorticity, h, 
becomes Infinite, F’ approaches F and G^ approaches G. Also as k tends to 
infinttv, C (kmh) approaches C(k). When h is zero. all the vortictty les in 
the plane of the airfoil, singulartties result, and the meaning of any obtained 
results. physically speaking, ts not clear [4] Loewv [4], and Couch [40] showed 
plots for F and G versus reduced frequency k. Loewy treats an infinite 
number of wakes and Couch explores a finite number including the special case 
of a single wake. These plots show the effects due to wake phasing, m , and 
wake spacing, h A case of special interest ts the wake shedding from one 
preceding blade only. lFor this case Loewy s lift deficiency. function, C^. reduces 


to the following lift function. C*, obtained by Wood and Couch [40] 


ወው SG: (4.6) 
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The propulsive force coefficient can be calculated and the results, obtained 5v 


Wood and Couch [40], using F* and G* are: 


For pure plunge : C, = wR AS (FP +(G")| 





For pure pitch. - C, 5 و‎ ን 
E 1 C 
MIC TT -a a 


D. EVALUATION OF THE TWO-FOIL PANEL CODE 


Several analyses are available in the aerodynamic !tterature “1hich permit a 


‘mited evaluation ot the two-foil panel code. These are: 
“round Effect on Airfoil Lift in Steady Flight 


5». . « ‘tka et al [44] analvzed the ground effect on a 2-D lal piate using 
the method of contormal transtormations. Figure + 2a shows the panet code 
results, Figure 42b reproduces Tomotika results. li seen that the two metnods 
are ın good agreement. Figure 4.3 a and Figure 4 2 b show the companson be- 
tween the result obtained bv Wegiev [46] tor a 19°% thick RAEIOI airroil in 
ground effect and the panel code prediction tor a NACA 0010 at zero angle of 
attack. The two results are again seen to be in good agreement considering the 
fact that the RAEIOI airfoil (whose coordinates are not known ) was approxi- 


mated by a NACA 0010 airfoil. 
2. Big....e Effect 


Glauert [45] defined a factor B which gives the lift reduction on the airtoil 
in an unstaggered biplane configuration compared to the lift on a monoplane at 
the same incidence. Table 4.1 gives the correction factor tor several non-dimen- 
sional biplane spacings. It is seen from this table and from Figure (4.4) that the 


two results are again in reasonable agreement. 
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Figure (4.2 b) Cl/Clo vs h/c ( From Ref. [44] ) 
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ligure (4.324) Local Velocity Distribution, V,/V vs A/C (Panel Code) 
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figure ( d. 3 b) Local Velocity Distribution, V,/V vs X/C (after Wegley [46]) 
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Figure (4.4) Correction Factor, B vs the Chord-gap Ratio h/c 
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TABLE 41 CORRECTION FACTOR FOR ^4 BIPLANE 
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0.50 | 0.770 0.730 
۸ 0 0.800 
| 1.00 | 0,902 0.855 | 
| 1.25 | 0.942 0.895 | 
150 | 0.97 0.920 | 
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Another ground interference study is due to Bagley [46]. As can be seen from 
the companson of the panel code (Figure 45a) with his predictions ( Figure 4.5 
bh), in ground effect, the lift is reduced over the upper surface and increased 
over the lower surface. As the airfoil approaches the ground, the flow between 


the airfoil and ground is reduced and hence the pressure is increased. 
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Cp vs X/C (SS.h/c2.5. AOA-11 deg.) 
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Figure (4.5 a) Pressure Distribution vs X/C ( Panel Code ) 
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E. ANALYSIS OF OSCILLATORY BLADE INTERFERENCE EFFECTS 


Having determined in Chapter IHI that pure rlapping ol a single nelicopter 
rotor blade 1s unlikely to produce a drag reduction sufficient to eXpiain the 
OH-6A HHC power reductions [1]. solelv. due t» tine Katzmavr cıfect, we 
investigate in this section the possibility ot favorable blade,wake interterence 
etfects. Consider the arrangement of the two airtoils snown in Figure (+6) The 
wake shed trom airtoil +2 impinges airfoil = 1. if the vertical spacing between 
the two airtoiis ralls below a certain \alue. depending on the +5 
oscillation ot airtoil #2. This case cannot be analyzea with the unsteady panel 
code because the vortices start to penetrate the second airo or because the 
wakes shed trom the two airtoils come in contact with each other. as shown 


in Figure (4 7). 


On the other hand. if the vertical non-dimensional spacing 1s two or more 
no problems occur. A second important parameter is the horizontal spacing be- 
tween the two airfoils. This parameter controls the phasing between the two 
wakes. Hence the unsteady two-foil code permits the anaivsis ot Loewy s ar- 
rangement for the special case of two wakes. This case is considered by Wood 
and Couch [40] in closed form and presented in Ref. [39]. Furthermore. the 
single airfoil case is recovered if the vertical spacing is chosen to be large 
because the effect of the first airfoil on the second must diminish with increas- 


ing distance between the two foils. 
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Figure (4.6) Wake Pattern (Airfoil #2 Plunging at High reduced Frequency 


and Close to Airfoil #1) 
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Figure (48). and Figure (4.9) show the computed litt and drag (thrust ) co- 
efficient tor airtoil #1 in the presence of airtoil = 2 Both airtoils are oscillating 
in plunge with an amplitude ot 0 14 and a reduced frequency ot 00617 The 
vertical spacing between the two airfoils is 200. In this case the lift and thrust 
coetficients are tound to be identical to the ones obtained from the single airfoil 
code. thus showing that a distance of 200 is sutfictent to recover single airfoil 


results. 


In Figures (4 10) to (4 21). on the other hand. the vertical spacing 1s reduced 
to 2. For these calculations. the time step ts 3 to.3 time units. The piunging 
amplitude is. again 0.14, the reduced frequency is 0.0617. The twelve figures 
4.10 to 421 show four different cases of phasing between the two wakes. In 
Figures (4.10) to (4.12), the wake pattern, lift and drag are shown for zero 
phasing, m = 0., Figures (4.13) to (4.15) show a phase of m=0.20833, Figures 
(4.16) to (4.18) show the case of m=0.25, and Figures (4.19) to (421) show 


the case of m=05. 


A closer inspection of these figures reveals the strong etfect of the wake 
phasing. For m =0, airfoil # 1 experiences a reduction in thrust compared to the 
single foil, whereas for m=0.5 airfoil #1 experiences a significant increase in 
thrust. These results are summarized in Figure (4.22) where the single airfoil 
results are also shown as well as the predictions by Wood and Couch [40] 
when specialized to the case of two interfering wakes only. The agreement be- 


tween the two predictions 1s encouraging. 
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Figure (4.14) 
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l'igure (4.16) Wake pattern ( plunge, m = 0.25 , h=2) 
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Wake pattern ( plunge, m=0.5 , h=2) 
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Similar results are obtained for the pitch case. The airfoils studied are again 
NACA 0007 foils. pitching about the quarter-chord point with an amplitude of 
| degree and a reduced frequency of 0.0617. Figures (4.23) to (4.32) display 
the same phasing values as in the plunge case. These figures and the summa- 
izing Figure (4.33) show again that values of m near zero and one produce 


increased drag compared to the single airfoil case, but values of in near 0.5 


produce thrust, whereas the single airfoil produces drag. Both analyses, the anal 
ysts by Wood and Couch [40] for the special case of two wakes only and the 
panel code, produce similar trends. However, the quantitative agreement is not 
as good as in the plunge case. The precise reason is not sufficiently understood 
at this time. It appears that the resolution of the suction peaks over the leading 
edge and therefore the computation of the thrust is more sensitive in the pitch 
case than in the plunge case. Furthermore, it must be remembered that Loewy's 
analysts for an infinite number of wakes and that by Wood and Couch [40] for 
two wakes are based on the flat-plate assumption, whereas the panel code requires 
a minimum blade thickness in order to obtain accurate results. A further evalu- 
ation of the panel code can be done for the drag produced by a single airfoil 

pitching about the leading edge. In this case the analytical theory of Garrick [5] 
and the computational results of Bosch [ 15] are available. It can be seen fron 
Figure ( 4.34 ) that the trends are again in agreement, but the code results deviate 
from Garrick’s and Bosch’s results with increasing frequency. These deviations 
must be attributed to the geometry differences ( flat plate versus NACA 0007 ) 
and the higher amplitude used in the panel code. Note that the pitching airfoil 


Starts to develop thrust only at a reduced frequency greater than about 1.8. 
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Figure (4.34) Average drag versus reduced frequency ( single airfoil ) 
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V. PROPULSIVE EFFICIENCY 
AND FLUTTER ANALYSIS 
A. INTRODUCTION 


In the previous two chapters it was shown that a single airfoil which 
pertorms harmonic plunge oscillations generates a certain amount of forward 
thrust (Katzmayr effect) and that this thrust may be enhanced by the oscillatory 
wake shed from the preceding blade if the frequency ratio or phasing m is in 
the range 0.2 < m > 0.7 approximately 0.5. It was also shown that a similar 
enhancement occurs for pitching blades. However, pitching airfoils develop thrust 
only for relatively high values of reduced frequency, k > 0.6 by Garrick 


results [5]. 


Furthermore. it is well known that a pitching airfoil may develop negative 
aerodynamic damping and thus experience single-degree-of-freedom flutter. As 
nown by Loewy [4], the possibility of flutter is enhanced by wake interference. 
This is an effect known as wake-induced flutter. It is therefore of interest to 
study the effect of various parameters on propulsive efficiency and on flutter 


instability in more detail. 
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B. PROPULSIVE EFFICIENCY 


The propulsive efficiency is defined as the ratio of the average work of 
propulsion to the work required to maintain the oscillation. If this factor 1s zero 
or negative. we expect energy to be extracted from the air and fed to the oscil- 
lating system. For a flat plate performing a flapping motion given bv h(t), where 
h(t) 2h, e'?' the result for the horizontal force per unit span. averaged over 


one cvcle. is given by [42]: 
D - -npbo^h" F^(k) -G'(k). (5.1) 


where F(k) and G(k) are the real and imaginary parts of the Theodorsen lift 
deficiency function C(k). For the case of pitching motion, the average horizontal 


force per unit span prs given by the following equation, (5.2): 


وی وب 
.—8a]—-—-—|--a‏ 
a AO s‏ 


k^ 


— 


1 1 Zu 
D= -npb'o a^ IF cd "Exi A +|5-F 


2 2 
k 





where a, is pitch amplitude and (a) is the position of axis of rotation. a= - | 
corresponds to the leading edge, a = +1 is the trailing edge. The propulsive 
efficiency then is given by : 
uet (5.3) 
: W 


W is the average work required to maintain the oscillation. For pitch, it is 


given by [42]: 


W = xpb'o Ua; (5 (5-a)-(2+5)]F JCEDE = (5p 
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For plunge, it is given by 
OS ae 
W = zpb'h o ui (k) = npbh^o UF (k) ay 


Ihen the efficiency factor, 1], can be expressed as, for plunge 


p 0 
O ٢ 


] > 
۳ F(k) 


and for pitch 


(5.7) 


Figure (5.1) shows the propulsive efficiency of a single plunging airfoil as 
a function of reduced fiequency. Slowly flapping airfoils reach very high values 
ol efficiency, approaching one. As the frequency of flapping is increased, the 
efficiency falls to 50 percent Pitching airfoils, on the other hand, start out with 
an efficiency of no more than SO percent at high frequency of oscillation which 
quickly falls to zero and to negative values with decreasing frequency. A negative 
propulsive efficiency means that no work is required to maintain the oscillation, 
but instead energy is extracted from the airstream. This implies a loss of damping 


and thus the possibility of single-degree-of-freedom flutter in pitch. 
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Figure (5.1) Propulsive Fficieney 1) versus l/k for a plunging airfon 
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Figure (5.2) shows the propulsive efficiency of an airfoil which pitches about 
the leading edge. The efficiency quickly drops to negative values as the fre- 
quency is reduced. As seen in Figure (5.2), infinitely large negative values are 
reached as the inverse frequency approaches 24, corresponding to a reduced fre- 
quency of 0.038. This is the flutter condition, as will be explained in more 


detail in the next section. 


Similar trends are found for other pitch axis locations, as seen in Figure 
5.3 ( pitch axis is at the quarter chord and mid-chord point ) The three quarter 
chord point is an exception. Here the propulsive etficiency 1s minus one at very 
high reduced frequencv, which falls to larger negative values with decreasing 
frequency ( Figure 5.4). However, it is interesting to note that a small shift of 
the pitch axis of only 0.5 % in either direction restores the efficiency to 0.5 


for large frequencies, from which it starts to drop with decreasing frequency. 


These results show that there is a fundamental difference between airfoil 
plunge and pitch motions. Plunge oscillations of a single airfoil always produce 
thrust. and the efficiency decreases as the frequency is increased. Pitch oscil- 
lations, on the other hand, produce thrust only at relatively high values of re- 
duced frequency, k > 0.6. This possibility of single-degree-of-freedom in pitch 


flutter is discussed in more detail in the next section. 
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ligure (5.2). Propulsive Efficiency 11 versus 1/k for an airfoil pitching 


about the leading edge 
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Figure (5.3) Propnlsive Efficiency versus 1/k for an airfoil pitching 


about the quarter - chord point and mid - chord point 
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Figure (5.4) Propulsive Efficiency n versus 1/k for an airfoil 


pitching about the three- quarter- chord point 
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C. SINGLE-DEGREE-OF-FREEDOM AIRFOIL FLUTTER 

Consider an airfoil hinged at its leading edge but elastically restrained from 
rotating about this axis by a torsion spring with constant Ką (ft-lb /rad). The 
airfoil is placed in a low-speed airstream so that the unstrained position of the 
spring corresponds to zero angle of attack qa, Figure (5.5). The equation of 


motion for this single-degree-of-freedom svstem is given bv : 
[a+K a@=M. (5.8) 
(L (ፈ \ 


M, is the aerodynamic moment due to a(t) and I, is the moment ot inertia 
about the leading edge. Assuming that  a(t)= a, e. then equation (5.8) can 


be written as 





— 
ER 
MS" | 


Da 
-|1- ES +m = 0. (5.9) 
tpb L a 


where o,, is the natural frequency of torsional vibration and is given bv : 


Ko 
a (5.10) 
T 


and m, is shorthand for the dimensionless aerodvnamic coefficient and can be 


written as [41]: 


^ 


d 


M 


19.1 7 1 
m, = 7 حم‎ ZENI -(L اھ : 6 ومو بسن‎ 
ነ د‎ 2 101 a a 9 2 h 2 
TPb 0 a e A 
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1 " 


where "a" is. the elastic axis position and at the leading edge a= -I m, 


[5 a complex number and it is a function only of reduced frequency, k — o b/U. 


[:quatton (5.9) can be split into real and imaginary. parts. 


| - | 
m | >. | 
Ref{m } = Re | MM eee 3) 
ya 0 (0 | 
apb 
E E e MK M T ee سر‎ (5412) 


Flutter occurs at that value of the reduced frequency where the imaginary 
( out-of-phase) part of the aerodynamic moment becomes zero. provided that the 


corresponding real (in-phase) part yields a non-tmaginary flutter frequency. 





Figure (5.5) An airfoil restrained to rotate about its leading edge 


in (two-dimensional flow. 


BUS 


Figure 5.6 shows the variation of the real and imaginary ( damping ) parts of 
the pitching moment for an airfoil pitching about the leading edge. It is seen 
that the pitch damping goes through zero at a reduced frequency of 0.038. This 
value agrees with the previously found value using propulsive efficiency (energy 


extraction ) considerations. 


Figure 5.7 displavs the companson between the pitch damping coefficients 
(note that the pitch damping coefficient is defined by Theodorsen [5] as: 

TS . E ١ | E | 

C = —— M. - E +a) O ር +a | L 


m 2 


The difference in Cm shown in Figure 4.8 is due to a difference between the 
Ma , L4, My and L, terms in. Equation 5.13 as determined by panel code 
( Riester [43] ) and Theodorsen's theory [3]. The imaginary part in Figure 5.7 
is obtained by Equation 5.11 computed by the panel code and Theodorsen's 
analysis. It is seen that the two computations are in reasonable agreement, 


recognizing the difference in geometry ( NACA 0007 versus flat plate ). 


Systematic variation of the pitch axis location shows that zero pitch damping 
ıs possible for axis locations upstream of the quarter -chord point. No flutter is 
possible. on the other hand, for axis locations at or downstream of the quarter 
chord point, Figure 5.8. Loewy [4] has shown that wake interference greatly 
increases the possibility of flutter. The unsteady panel code provides an oppor- 
tunity to investigate and to extend this well-known finding. This is a problem 


which should be studied in more detail. 
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l'igure (5.6) Variation with reduced freqnency k of the real and imaginary 
parts of the non-dimensional acrodynamic moment m, due to pitching of 


an airfoil about its leading edge in an compressible flow. 
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l'igure (5.7) Comparison of Cm (Imag.) vs Kp for NACA0007 Oscillating 


m pitch abont the leading edge for o =1° and N =200 panels. 
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l'igure (5.8) Imag. (my ) vs k (for an airfoil pitching at 
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VI. RESULTS AND CONCLUSIONS 
A DISCUSSION OF RESULTS 


This research was undertaken. to study the effect of unsteady aerodynamics 
on reductions of main rotor shaft torque and engine power measured when 
the HHC system was applied to the OH-6A rotor. In Chapter III we study pro- 
pulsive force resulting from pure plunge oscillation of an airfoil in inviscid, 
incompressible flow. The results of flat plate theorv, Theodorsen [3] and Garrick 
[5], are used to validate panel code results. Comparison with Garrick showed 
good agreement in pure plunge, especially, for lower values of reduced frequency 


and amplitudes of oscillation. 


Figures 3.5 to 3, show the average drag coefficient versus plunging 
amplitude for three values of reduced frequency, k= 0.1, 1. and 2 respectively. 
It is seen from Figure 3.5. ( reduced frequency, k is 0.1) that agreement 
between the panel code and Garrick [5] is very good for a range of non- 
dimensional plunging amplitude up to 15 to 20 % of the blade chord. On the 
other hand, for higher values of k, the two results are in good agreement for 
a range of plunging amplitude up to approximately 10%. Single airfoil code, was 
used to compute the time history of the drag (thrust) coefficient of the NACA 


0015 airfoil ( the blade section of the OH-64A ). 
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Figure 3.12 and Table 5.3 show the average drag coefficient Versus reduced 
frequency for an NACA 0015 airfoil which executes plunging oscillation with 
|" plunging amplitude. The OH-6A data presented in Table 3.1. Wood [42], 
and Figures 3.8 to Figure 3.11 from Ref. [27], were used to estimate the 
propulsive force per blade obtained due to plunging oscillation tor | P frequency 
case. In these figures. the difference between the deflections when the HHC 
is "on" and " off " represents the plunging amplitude. The controller phase is 
also significant. Phase angle of 90° or 120° are those where helicopter vibra- 


tlons were most sever. whereas at phase angle of 300° or 330". the helicopter 


vibrations were significantly reduced. 


We look for the greatest power benetit to occur at the phase angle of 
lowest vibration, that is 300° or 330°. Table 3.2 shows the amount of propulsive 
force and torque generated at each blade segment. The results of this table 
showed that the contribution of 1P (first harmonic) were 3 horsepower. Wood 
et al [27], applied the Garrick equation and considered the effect of blade plunge 
response in relation to the second to twelfth harmonic loads. He found for 2 P 
harmonic only ( HHC "on", 330° phase at 60 knots ), the power gained was 
about 11.3 horse power. For more detail. the reader is encouraged to study Table 


«የ ክብ. 


I Wood E. R., Hıgman. J. and Ramesh Kolar. " Higher Harmonic Control 
Promises Improved Dynamic Interface Operations" AGARD Proceedings of 
the 78 [ከ Flight Mechanics Panel, Mav 1991, Seville. Spain. 
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Chapter IV, considers the influence of the rotor wake interaction and examines 
whereas such interaction might cause favorable effects. First, the case of an iso- 
lated airfoil ( single airfoil oscillating in plunge) with absence of wake was 
considered ( or the wake located at great distance from rotor ). This is the "Katz- 
mayr effect" case where the resulting drag force is a propulsive force. It was 


the primary case considered by Garrick [5]. 


The advantage of the Panel Code is that it provides us with time histories 
of unsteady lift. drag, and the wake trajectorv bv which we can better understand 
and analyze the problem. n each case. the wake pattern was plotted in order 
to be sure that the wake produced by the second airfoil would not interfere with 
the reference airfoil. Also the lift and drag(or the propulsive force) were plotted 


versus non-dimensional time, T* ( T* = t.U/c ). 


The closed form solutions are based upon a very thin flat plate airfoil 
undergoing very small motions when we model it in the unsteady panel code. 
For the calculations an NACA0007 airfoil was selected. The airfoil sets at zero 


mean angle of attack and the reduced frequency, k. 1s = 0.0617. 


For the OH-6A, the value of reduced frequency, k at 70 % radial station, 
was calculated and found to be .1234 for the 32 Hz exciting frequency. The 
reference length used by the panel code is the blade chord or 2b, so the 0.1234 
reduced frequency corresponds to 0.0617 per the panel code definition. Non- 
dimensional wake spacing, for single airfoil analysis, was h=200 ( mathemat- 
ically equivalent to infinite spacing, or we say that the interaction between the 


two airfoil is very small). 
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The amplitude of plunge oscillation was h,= 0.14 ( bv the panel code 
.07, for the same reason ). The panel code results are compared to other refer- 
ences in order to know the code limits. The comparison of the code results to 
that obtained by Ref. [44], is shown in Figure (42). This figure shows the 
relative lift coetficient of an airfoil (a wing of relatively long span) as it ap- 
proaches the ground for different values of angles of attack. The tigure. showed 
that the results estimated bv the panel code and the results obtained by Ref. 
[44], are very close. For h/c < 0.5 theresults obtained bv the panel code 


are suspected. 


Another comparison is shown in Figure (4.4). This figure and Table 4.1 
show the correction factor, B, versus the chord-gap ratio in a biplane contigura- 
tion. The correction factor, as defined by Ref. [45], is the reduction of the lift 
coefficient of an unstaggered biplane compared with that of monoplane at the 
same angle of incidence. On Figure (4.3) and Figure (4.5), the local speed ratio 
and pressure coefficient vs X/C are plotted and showed very good agreement 


with Wegley [46] and Bagley (1959), respectively. 


The panel code is suspect whenever the wake generated by an airfoil 
impinges on another airfoil. The numerical panel code method also encounters 
difficulty with discontinuities. For these reasons it is important to carefully check 


the time history plots after each computer run. 
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In Figure (4.8), the time history of the lift coefficient. Cj. vanation is 
plotted. Shown in Figure (4.9 is the drag coefficient, Cy versus the non- 
dimensional time T*. We observe from Figure (4.8) that the lift varies as the 
frequency of oscillation, with a mean lift value of zero as we expect. Amplitude 
of oscillatory lift was found to be+ 0.0510. Figure (4.9) clearly shows the " 
Katzmayr effect". We noted that the calculated drag varies at twice the plunge 
oscillation frequency of the airfoil. Also, observe that the airfoil is generating a 
propulsive force as indicated by the mean value of drag coefficient, which is 


seen to have a negative value. 


Looking to the computer printout at, at the non-dimensional time T* = 
0.0, steady state, the drag generated by the first airfoil was found to be 
Cas, = -000131 for NACA0007, which is the airfoil type that was used for 
most of our analyses at zero angle of attack, and a total number of panels 
N= 100 ( 1.e 50 panels at the upper surface and 50 panels at the lower surface). 
For NACA0009, (this airfoil was used to compare the code results with that of 
Bosch's for a pitching case), the steady state drag C4,, was found to be .000196 
at N — 100 and .000113 at number of panels N =150. 

It is known that, for an airfoil at zero angle of attack, in an inviscid 
incompressible flow field, steady state drag is zero, so in our calculations that 
steady state drag was subtracted from the drag generated at each time step. Shown 
in Figure (6.1) is the steady state drag coefficient vs angle of attack [37]. Figure 
(6.18) shows the variation of steady state drag coefficient vs angle of attack for 
three different airfoils; NACA0003, NACA0007, and NACA0012 at total number 
of panel of 200. 
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In this graph, it is shown that. for the thin airtoii NACAOOOS at zero angle 
of attack. the drag coetficient is very small but as the angle of attack increases. 
the drag coetficient highly increases more than in the cases of NACA0007 and 
NACA0012. Figure (6.1b) shows the etfect of increasing the total number of pan- 
els N on the drag coetficient of the airfoil NACAO0012. As the total number 


of panels increases. the steadv state drag coefficient accumulated error decreases. 


Figures (4.10) through Figure (4.21), show the same case as above. that 
of an airfoil oscillating in plunge, but with the introduction or a wake layer 
in the near proximity of the reference airfoil (non-dimensionai wake spacing, h 
=2 ) for wake phasing, m = 0, 0.20833, 0.25, and m = 0.5. These correspond 


to phase angles of 0% , 75%, 90% and 180° , respectively. 


We see that the phase angles relate to the phase relationship between the 
reference airfoil and the wake layer of shed vorticitv. These phase angles were 
obtained in the panel code by establishing the proper time phase between the 
Initiation of the wake vorticity and the oscillation of the reference airfoil. This 
time phase is in turn a function of the frequency of oscillation which is directly 
related reduced frequency, k. For the results presented here, for k —.0617, this 
corresponds to 50.92 time units for m=0.0 (0° phase ); 61.53 time units for m 


= 0.2083 (75" phase); and 76.38 time units for m — 0.50 ( 180? phase ), etc. 
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If the reduced frequency. Kk, changes. this corresponds to a change in 
frequency (at given chord length and free stream velocitv ). For k =0 1234. this 
correspond 50.92 time units for m = 0.0, 56.23 time units for m — 0.2083, and 
63.65 time units for m — 0.5, etc. 

Figure (4.10) shows the wake position for m= 0.0, the case where the 
wake shed from the second airfoil at h=2 is in phase with wake shed from 
the reference airfoil itself. In Figure (4.21), the time history of lift variation 
is shown. 

This 1s indicated by the solid line for Cl, . The dashed line. Cl», indicates 
the variation of lift on asecond airfoil located at the point of initiation of the 
wake layer at h=2. Our discussion will focus airfoil, , indicated by the solid 
line in the figures mentioned above. For an airfoil oscillating about zero angle 
of attack, the mean lift is zero. We found that, while the mean lift is zero, the 
amplitude of lift oscillation is smaller in the presence of the wake at zero 
phasing. From these figures. we conclude that the effect of wake vorticity at 
zero phasing is to reduce the " Katzmayr effect". 

Comparison of oscillatory lift values of the 75? wake phasing case, with 
the case of no wake or single airfoil, shows that they are identical. since the 
values at 75 degrees phase are the same as that for infinite wake. Table 6.1 


proves this result. 
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The time history of drag is seen to varv as twice the plunge oscillation 
frequency of the airfoil. We also observe that the airfoil is generating propulsive 
force as indicated by the mean value of drag, which is seen to have a negative 
value. Results of 75 degree wake phasing case, contirm that there is a phase 
angle. m, at which the effect of the wake vorticity on the reference airfoil gives 
the same results as when the wake is removed . It has been shown at a phase 
angle of zero degrees that the shed wake vorticity increases the drag on the 


oscillating rete:ence airfoil. i.e. diminishes the "Katzmavr effect". 


Conversely we found that at a phase angle of about 180 degrees. there is a 
significant decrease in drag at the reference airfoil due to the wake vorticity. 
This implies an intermediate phase angle at which the drag is identical to the 
"Katzmayr effect". It turns out that such a phase angle exists. In fact, at this 
phase angle the lift, drag and pitching moment are found to be identical to the 
case of no wake, and match the "no wake" values identically. 

Table (6.1) show the aerodynamic coefficients CI(1), CI(2), Cm(1), Cm(2), 
Cd(1), and Cd(2) vs. Non-dimensional time. This table is a part of one of the 
output files. It shows these aerodynamic coefficient vs time from the nondimen- 
sional time unit 113.3 until 159.1. The phase angle at which this occurs is 75 


degree or m — 0.20833, for the case considered here. 
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Table (6.1) Aerodynamic coefficients vs Non-dimensional time (from time 
step 113 to 159) for the case of wake phasing, m=75 degree, wake spacing 


l= 2. reduced frequency k = 0.0617 and plunging amplitude hbar = 0.14 ) 


TINE eno) CL(2) CM(1) CM(2) CD(1) CD(2) 
113.290932 0.006952 0.014965 -0.001387 -0.003432 -0.000011 -0.000039 
114.563866 -0.001073 0.007130 0.000665 -0.001427 0.000000 -0.000011 
115.836800 -0.009081 -0.000881 0.002701 0.000614 -0.000009 0.000000 
17.109734 -0.016877 -0.008871 0.004674 0.002639 -0.000036 -0.000008 
11A.3A2668 -0.024266 -0.016644 0.006534 0.004599 -0.000079 -0.000035 
119.655602 -0.031067 -0.024005 0.008236 0.006446 -0.000133 -0.000078 
۱20.928515 -0.037111 -0.030777 0.009738 0.008134 -0.000194 -0.000132 
122.201469 -0.042747 -0.036791 0.011001 0.009621 -0.000255 -0.000193 
123.474403 -0.046347 -0.041898 0.011994 0.010872 -0.000310 -0.000254 
124.747337 -0.049307 -0.045974 0.012693 0.011855 -0.000355 -0.000309 
126.020271 -0.051051 -0.048918 0.013077 0.012546 -0.000384 -0.000353 
127.293205 -0.051534 -0.050659 0.013140 0.012929 -0.000394 -0.000382 
128.566132 -0.050746 -0.051152 0.012878 0.012994 -0.000386 -0.000393 
129.839066 -0.048701 -0.050385 0.012296 0.012738 -0.000358 -0.000384 
131.112000 -0.045454 -0.048379 0.011412 0.012169 -0.000315 -0.000358 
132.384933 -0.041083 -0.045182 0.010245 0.011301 -0.000260 -0.000315 
133.657867 -0.035698 -0.040872 0.008825 0.010155 -0.000199 -0.000261 
134.930801 -0.029433 -0.035556 0.007188 0.008758 -0.000138 -0.000200 
136.203735 -0.022446 -0.029364 0.005375 0.007146 -0.000082 -0.000139 
137.476669 -0.014910 -0.022450 0.003431 0.005358 -0.000038 -0.000084 
138.749603 -0.007014 -0.014982 0.001404 0.003437 -0.000009 -0.000040 
140.022537 0.001046 -0.007147 -0.000655 0.001432 0.000001 -0.000011 
141.295471 0.009072 0.000865 -0.002696 -0.000608 -0.000007 0.000000 
142.568405 0.016865 0.008856 -0.004668 -0.002634 -0.000035 -0.000008 
143.841339 0.024230 0.016627 -0.006521 -0.004593 -0.000078 -0.000035 
145.114273 0.030990 0.023991 -0.008212 -0.006441 -0.000133 -0.000078 
)46.387207 0.036981 0.030761] -0.009700 -0.008128 -0.000194 -0.000132 
147.660141 0.042057 0.036775 -0.010948 -0.009616 -0.000255 -0.000193 
148.933075 0.046093 0.041883 -0.011926 -0.010867 -0.000310 -0.000254 
150.206009 0.048996 0.045960 -0.012610 -0.011850 -0.000354 -0.000309 
151.478943 0.050693 0.048903 -0.012985 -0.012541 -0.000382 -0.000353 
152.751877 0.051145 0.050643 -0.013039 -0.012923 -0.000392 -0.000382 
154.024811 0.050343 0.051136 -0.012775 -0.012988 -0.000384 -0.000393 
155.297745 0.048306 0.050371 -0.012196 -0.012733 -0.000357 -0.000384 
156.570679 0.045084 0.048364 -0.011319 -0.012164 -0.000314 -0.000357 
157.843613 0.040757 0.045167 -0.010163 -0.011296 -0.000260 -0.000315 
159.116547 0.035427 0.040858 -0.008757 -0.010150 -0.000199 -0.000260 
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The proper phasing is achieved by positioning the reference airfoil 76.38 
time units to the rig. of the second airfoil, the airfoil that generates the wake. 
We found in the case of the plunging airfoil that the wake at zero degree phase 
Increases the drag force. while the wake at 75 degree phase gives us close 
results to "no wake at all" case. So, we expect that the wake at 180 degree to 
increase the propulsive force. Shown in Figure (4.22), is the variation of the 
average propulsive force coefficient versus wake phasing, m. The circles in the 
graph represent panel code results for the parameters mentioned above. 

Superimposed on the plot of the figure are the analytical closed form 
results of the Loewy theory as modified for one wake, as shown in Figure 
(4.22) from Wood et al [38]. Note the good agreement between the two sets 
of results. It is likely that the agreement can be further improved by increasing 
the number of panel from 100, used for these calculations, to 200 or 200 more. 

As was done previously for an airfoil oscillating in plunge, we will now 
consider the case of an airfoil oscillating in pitch with particular emphasis on 
drag or propulsive force. We will first look at the case without presence of 
layers of adjacent wake, then look at the wake effect including the important 
consideration of phasing. We will begin by considering the numerical results and 
comparing these results with Garrick, Bosch, for "no wake" case, and then com- 


pare with Loewy’s theory when the wakes are included. 
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Considered is the case of an NACAO007 airtoil set at zero degree angle 
of attack. The airfoil is oscillating with an amplitude | degree. The reduced 
frequency is taken at k=0.0617 and the non-dimensional wake spacing was 
h =200, which is numerically equivalent to having the wake at infinitv. The 6 
history of the lift variation, Cl, and the time history of drag variation. Cd, reveal 
the following. We observe that the lift varies as the frequency of oscillation. 
with a lift mean value of zero as we expect. 

In contrast to the lift time history, the drag time history indicates that the 
calculated drag varies at twice the frequency of the pitch oscillation. The am- 
plitude is given by Cd= 0.000192. The mean drag value 15 found to be Cd= 
0.00009 and the mean lift is of course 0.0 at zero angle of attack. 

For an airfoil oscillating in pitch with the wake layer near the reference 
airfoil at non-dimensional wake spacing, h=2. The case of zero degree phase, 
where the wake below the airfoil at h=2 is in phase the wake shed from the 
reference airfoil itself. Review of the results indicates a reduction in the ampli- 
tude of lift variation for the m=0.0 case when compared to the case with no 
wake at all (wake at infinity ). We see that the lift amplitude decreases. We also 
see an increase in mean steady drag accompanies the reduction in lift amplitude, 
for the case of zero wake spacing (m=0). Proper phasing is achieved at m 
=0.5. At this value, we found that the lift enhancement in this case of pure 


pitch, is similar to what we found in the plunge case. 
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B CONCLUSIONS 


The most important conclusions that can be drawn from this research are; 
1. For drag, in both cases of pure plunge and pure pitch, wake phasing can 
increase or decrease the steady component of drag acting on the airfoil. The 


increase in the steady component of propulsive force due to wake phasing in the 


case of pure plunge was sufficient to significantly enhance the " Katzmavr 


effect value of propulsive force. 


2. For an airfoil oscillating in plunge or pitch, tor optimum reduction in drag 
and enhancement in vibratory lift, the phase angle of the wake vorticitv in the 
single wake case with respect to motion of the reference airfoil should be about 


180 °. 


3. For an airfoil oscillating in plunge or pitch, the largest value of steady drag 
and smallest value of oscillatory lift occurs when the phase angle of the wake 


vorticity with respect to motion of the reference airfoil is about 0°. 


4. For lift, in both plunge and pitch in the single wake case. wake phasing can 
increase or reduce the oscillatory lift acting on the airfoil. Where in the absence 
of wake vorticity, the effect of oscillations decreases the lift (lift deficiency ). 


With the wake present we also observe lift enhancement ( lift efficiency ). 


5. There exists a phase angle, 75° for the case, (reduced k = .0617, wake 
spacing, h = 2), where the effect of wake vorticity on the reference airfoil is 
identical to the case of the wake at infinity ( represented in the code by wake 


h= 200), or no wake vorticity at all. 
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Applying these findings to the OH-6A case, we can add the tollowing 


conclusions : 


n 


6. Measured reductions in power are feasible when ıhe "Katzmayr" effect and 


thesis - reported wake enhancements are included. 


7. Measured higher harmonic control results show least benefit at 100 knots 


due to the wake being transported greatest distance from rotor at this speed. 


8. Poor repeatability of measured open loop performance data is to be expected 


due to shifts in wake position `. 


9. For helicopters, it is standard rotor blade design to locate the center of 
gravity, aerodynamic center, and elastic axis at the quarter chord point to avoid 
flutter. This was illustrated in the analysis. Results of this thesis showed that 
flutter is impossible for this condition without wake interference. However, 
Loewy’s results introduce the possibility of wake induced flutter. Future investi- 


gations, therefore, will have to explore this possibility in more detail. 





LER. Wood. Max F. Platzer. Ahmed Abourahma, Mark A. Couch 


" On the Unsteady Aerodynamics of Higher Harmonic Control. " 
Paper No. C 17, Nineteenth European Rotorcraft Forum, Ccmobbio 
(Como), Italy, September 1993. 
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C RECOMMENDATIONS 


A new flight test program. directed at studying helicopter performance, 1s 


needed. The NASA-Army OH-6A program was dedicated to vibration testing. 


Continued full scale HHC tests in the NASA Ames 40° x 80° wind tunnel is 
encouraged. In addition. it 1s suggested that tests be initiated in small university 


wind tunnels to verify the findings of this thesis. 


Controlled whirl tower testing 1s also recommended. since this allows careful 


measurement and control of performance parameters. 


A wind tunnel program is needed to explore similar cases and to verify the 


present research findings. 


. The panel method is a significant analysis tool that can be applied to study and 
solve similar problems. It provides an excellent method for studying single or 


two-degree-of-freedom flutter problems. 
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